HD-817*  672  0VN88IC  HMD  THERMODrNRNIC  CHHRRCTERISTICS  OF  8 

HICROBURST-PRODUCINQ  STORM.  .  <U>  SIR  FORCE  INST  OF  TECH 
HRIQHT-P8TTERSON  NFS  OH  R  0  HUGHES  1986 


UNCLASSIFIED  8FIT/CI/NR-86-69D 


I 

MICROCOPY  RESOLUTION  TEST  CHART  * 

NATIONAL  BURIAU  Of  STANDARD^  !<*,«  i 


AD-A170  673 


Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
FJF.FORE  COMPLETING  FORM 


7.  AUTHORfJj  , 

Robert  G.  Hughes 


performing  organization  name  and  address 

AFIT  STUDENT  AT:  Saint  Louis  .University 


8.  CONTRACT  OR  GRANT  NUMBERfsJ 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 


CONTROLLING  OFFICE  NAME  ANU  ADDRESS 

AFIT/NR 

WPAFB  OH  45433-6583  I  is.  number  of  pages 

214 _ 

14.  MONITORING  AGENCY  NAME  S  ADDRESS(iY  dittcrent  Ir oin  Controlling  Oilier )  I  15.  SECURITY  CLASS,  (o I  this  report) 


12.  REPORT  DATE 

1986 

13.  NUMBER  OF  PAGES 


UNCLAS 


15  a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


DISTRIBUTION  ST  ATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 


18.  supplementary  NOTES 

APPROVED  FOR  PUBLIC  RELEASE:  IAW  AFR  190-1 


IS.  KEY  WORDS  (Continue  on  reverse  side  if  ncce s  sary  and  identify  by  block  number) 


DTIC 

jZLECTE 

AUG  1  2  1981 

l^Ir 


ETON  E.  WOLAVER  &  4l/6  f% 
Dean  for  Research  and 

Professional  Development 
AFIT/NR 


20  ABSTRACT  (Continue  on  rovers *t  side  if  neces sary  and  identify  b>*  block  number) 


ATTACHED. 


86 


1  JAN  73  1473  EDITION  OF  I  NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flVirn  Dara  Enlrrmd) 


Z  W: 


Some  dynamic  and  thermodynamic  structures  of  a 
microburst-producing  storm  which  occurred  on  14  July 
1982  in  Colorado  were  studied  in  detail.  Dual-Doppler 
radar  data  collected  during  the  Joint  Airport  Weather 
Studies  (JAWS)  project  at  Denver's  Stapleton 
International  Airport,  were  objectively  analyzed  to 
produce  a  3D  wind  field.  The  analyzed  domain  had 
dimensions  of  10  km  by  10  km  by  8.5  km  centered  on  the 
microburst.  Vertical  velocities  were  computed  by 
integrating  the  anelastic  continuity  equation  downward 
from  the  storm's  top.  A  variational  approach  was  then 
employed  to  adjust  the  derived  three-dimensional  wind 
components.  Subsequently,  fields  of  deviation  pertur¬ 
bation  pressure,  density  and  virtual  temperature  were 
recovered  from  a  detailed  wind  field  using  the  three 
momentum  equations.  These  retrieved  fields  were  sub¬ 
jected  to  internal  consistency  checks  to  determine  the 
level  of  confidence  before  interpretation.  '■f-- 

Our  study  demonstrates  that  thermodynamic  re¬ 
trieval  is  feasible  even  when  the  storm  intensity  is 
non-severe.  Variational  adjustment  substantially 
reduces  errors  in  veitical  velocity  fields.  Results 
show  that  the  microburst  being  investigated  occurs 
during  the  decaying  stage  of  a  storm.  It  is  embedded 


2 


within  the  high  reflectivity  region  with  heavy  preci¬ 
pitation.  A  strong  downdraft  impinges  upon  the 
surface  producing  a  stagnation  mesohigh  inside  the 
microburst.  This  high  is  accompanied  by  low  pressure 
in  the  strongest  outflow  regions  forming  a  pronounced 
horizontal  perturbation  pressure  gradient  outward  from 
the  center.  Such  pressure  patterns  are  in  good  agree¬ 
ment  with  the  surface  observations  in  similar  cases. 
The  outflow  region  extends  from  the  surface  to 

approximately  1  km  height  AGL  with  maximum  divergence 

1 

in  excess  of  10  s  .  In  the  middle  troposphere,  high 
pressure  forms  in  the  upshear  side  of  the  main  updraft 
with  low  pressure  on  the  downshear  side  due  to  dynamic 
interactions  between  the  updraft  and  the  sheared 
environmental  wind.  Favorable  perturbation  pressure 
gradients,  evaporation  and  precipitation  loading  are 
responsible  for  initiating  and  sustaining  a  downdraft 
which  produced  the  microburst  at  low  levels. 
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curve  B  variational  analysis  with 
functional  J  ;  Curve  C  filtering  using 
the  Cressman  weighting  function.  The 
response  curve  D'  obtained  from  Barnes 
(1973)  scheme  with  V  =0.3  and  R=3.6  km  is 
denoted  by  the  symbol  x. 

Fig.  A8  Comparison  of  response  curves  obtained 
from  Barnes  (1973)  scheme  (D')  with 
V=0.3  and  R=2.5  km  and  theoretical 
formulas  derived  by  Testud  and  Chong 
(1983)  for  n-1  (Gl)  and  n=2  (G2). 

Fig.  B1  Relationship  between  "true"  functional 
solution  and  neighboring  function  for 
method  of  calculus  of  variation. 


CHAPTER  1 


Introduction 

In  early  studies  of  the  Thunderstorm  Project, 
Byers  and  Braham  (1948)  pointed  out  the  co-existence 
of  strong  updrafts  and  downdrafts  in  thunderstorms 
during  the  mature  stage  with  vertical  velocities 
reaching  40  ms.  They  suggested  three  mechanisms 
important  in  the  development,  enhancement  and  mainte¬ 
nance  of  the  strong  downdraft;  1)  entrainment  of  out¬ 
side  air  into  the  rear  of  the  storm,  2)  evaporation 
of  precipitation  into  this  drier  entrained  air,  and 
3)  aerodynamic  drag  effects  of  precipitation.  For  a 
severe  thunderstorm,  Newton  (1963)  added  to  these 
mechanisms  the  effect  of  forced  lifting  created  by 
this  downrush  of  colder  air  spreading  out  under  the 
storm.  This  forced  lifting  provides  a  continuing 
energy  source  in  the  form  of  warmer  moist  inflow. 
Both  Byers  and  Braham  (1948)  and  Newton  (1963)  along 
with  many  other  investigators  recognized  the  need  to 
understand  this  downdraft  and  its  role  in  the  thunder¬ 
storm  life  cycle.  Byers  and  Braham  (1948)  stated  in 
their  conclusions  that  the  distribution  of  temperature 
and  hydrometeors  within  the  thunderstorm  remained  to 
be  studied.  Newton  (1963)  in  his  summary  stated  that 


the  Thunderstorm  Project  of  Byers  and  Braham  had  pro¬ 
vided  many  new  hypotheses  which  needed  to  be  verified. 

Studies  of  the  internal  structure  of  the  thunder¬ 
storm  using  the  radar  as  a  remote  sensing  device  has 
been  in  progress  since  the  operators  of  the  Illinois 
State  Water  Survey  radar  spotted  the  first  hook  echo 
in  1947.  Browning  and  Ludlam  (1962)  documented  the 
internal  structure  and  kinematic  features  of  the 
Wokingham  supercell  hail  storm.  However,  they  were 
unable  to  detail  dynamic  and  thermodynamic  information 
due  to  the  restrictions  of  the  radar.  Others 
attempted  alternative  techniques  to  gain  insight  into 
the  internal  structure  of  storms.  Atlas  and  Wexler 
(1965)  suggested  a  modification  to  the  beam  pattern  of 
a  conventional  radar  which  would  allow  horizontal  wind 
fields  to  be  extracted.  Davis-Jones  and  Henderson 
(1973)  launched  in-storm  rawinsondes  and  Sinclair 
(1978)  attempted  aircraft  penetration.  All  these 
attempts  proved  to  be  too  expensive,  unsafe,  or 
unreliable  in  the  continuing  investigation  of  the 
thunderstorm. 

The  advancement  of  Doppler  radar  techniques  in 
deducing  air  motion  within  convection  storms  has  added 
the  needed  tool  necessary  to  increasing  our  under¬ 
standing  of  the  kinematic  and  dynamic  structure  of  the 
storm.  Multiple  Doppler  radars  can  provide  three 


dimensional  winds  inside  the  storm;  for  example,  see 
studies  by  Fille',  Justio  and  Rogers  (1963)  and  Armijo 


(1969).  Lhermitte  (1970)  was  the  first  to  use  the 
Doppler  radar  to  observe  horizontal  wind  fields  within 
a  convective  storm.  Ray  e_t  a_l.  (  1975)  were  the  first 
to  study  a  severe  storm.  Their  results  of  low  level 
wind  features  essentially  verified  features  described 
by  Ludlam  (1963)  and  Newton  (1963).  Since  the 
pioneering  work  of  Ray  e_t_  £l_.  (1975),  many  other 
researchers  have  used  the  multiple  Doppler  radar  as 
the  primary  tool  in  the  investigation  of  kinematic 
features  of  severe  thunderstorms,  e.g.,  Kropfli  and 
Miller  (1976),  3randes  (1977a,  1977b,  1973),  etc. 

Dual-Doppler  radar  measurements  cannot,  however, 
provide  a  direct  measure  of  thermodynamic  and  pressure 
information.  Fluctuations  in  pressure  and  temperature 
(potential  temperature)  are  important  in  understanding 
the  structure  and  internal  dynamics  of  a  convective 
storm  These  fluctuations  are  only  a  few  degrees  in 
temperature  and  a  few  millibars  in  pressure  and  cannot 
accurately  be  measured  by  conventional  upper  air 
sounding  techniques.  Some  type  of  indirect  measure¬ 
ment  method  is  necessary  to  recover  these  thermo¬ 
dynamic  fluctuations  from  measurable  fields,  such  as  a 


detailed  wind  field. 


Recently  Gal-Chen  (1973)  proposed  a  retrieval 
technique  in  which  the  momentum  equations  and  the  ob¬ 
served  three-dimensional  wind  field  are  employed  to 
determine  the  density,  pressure  and  temperature  fluc¬ 
tuations.  This  method  has  been  employed  on  model 
generated  data  and  observed  wind  fields  by  Gal-Chen 
(1978),  Hane  e_t  a_l .  (1981),  Lin  and  Pasken  (1982a, b), 
Pasken  and  Lin  (1982),  Brandes  0934),  Hane  and  Ray 
(  1985),  and  Lin  e_t  al .  (  1985). 

Since  the  method  proposed  by  Gal-Chen  (1973)  is 
indirect,  its  effectiveness  to  a  large  extent  depends 
on  the  degree  of  accuracy  in  the  input  data  (Hane  et 
al.,  1931).  Thus,  verification  of  the  recovered 
fields  of  pressure,  density  and  temperature  fluctua¬ 
tions  is  critical.  While  direct  verification  is 
desirable,  they  are  difficult  to  obtain  especially 
within  the  storm.  Consequently,  some  indirect  verifi¬ 
cation  is  needed  in  order  to  validate  the  method  and 
the  recovered  fields.  u  .1-Chen  and  Hane  (1981)  sug¬ 
gested  three  indirect  verification  techniques:  1) 
numerical  simulations,  2)  momentum  checking,  and  3) 
time  continuity.  Hane  e_t  a_l_.  (1981)  applied  this 
retrieval  technique  to  numerical  simulations  as  a 
guide  to  verification.  Momentum  checking  was  employed 
in  studies  by  Gal-Chen  and  Kropfli  (1984),  Hane  and 
Ray  (  1985)  and  Lin  e_t  a  1 .  (  1985).  Time  continuity 


checking  was  reported  in  a  study  by  Gal-Chen  and 
Kropfli  (1984). 

The  aforementioned  studies  showed  that  the 
retrieval  method  is  feasible  for  recovering  perturba¬ 
tion  pressure,  virtual  temperature  or,  virtual  poten¬ 
tial  temperature  within  a  convective  storm.  The 
degree  of  accuracy  of  the  observed  winds  appears  to  be 
the  k^y  factor  in  determining  the  effectiveness  and 
validity  of  the  method.  As  Gal-Chen  and  Hane  (1981) 
concluded  "Studies  using  real  data  are  encouraging. 
It  appears  that  pressure  fluctuations  within  30” 
accuracy  can  be  retrieved.  Retrieval  temperature  and 
pressure  fields  demonstrate  the  importance  of  non¬ 
hydrostatic  pressure  gradients  in  convective  dynamics 
and  appear  to  display  continuity  and  coherence.  Firm 
verification  of  these  results,  as  well  as  further 
technique  developments  (particularly  in  retrieval  of 
vertical  velocities),  should  be  pursued  with  vigor. 
Almost  everything  in  this  vast  and  difficult  problem 
still  remains  to  be  done." 

The  Joint  Airport  Weather  Studies  (JAWS)  project 
was  jointly  conducted  by  scientists  from  the  National 
Center  for  Atmospheric  Research  (NCAR)  and  the 
University  of  Chicago  in  the  summer  of  1982  at 
Denver's  Stapleton  International  Airport.  Primary 
objectives  of  the  project  include:  1)  the  research  on 
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fine-scale  structure  of  thunderstorm  dynamics  and 
kinematics  in  the  vicinity  of  the  airport,  2)  the 
effect  of  thunderstorm  low-level  wind  shear  on  air¬ 
craft  performance,  and  3)  the  development  of  real-time 
testing  of  low-level  wind  shear  interaction  and 
warning  techniques  and  displays,  see  McCarthy  £_t  al . 

(  1932)  . 

Fujita  (1978)  defined  the  term  "microburst"  as  a 
small  downburst,  less  than  4  km  in  outflow  diameter, 
with  peak  winds  lasting  only  2  to  5  minutes.  A  micro¬ 
burst  could  occur  in  areas  of  precipitation  or  no 
precipitation.  It  induces  a  sudden  outflow  of  damaging 
horizontal  winds  near  the  ground.  Fujita  and  Vakimoto 
(19S3)  further  noted  that  this  outburst  may  be  sensed 
at  the  surface  within  a  few  minutes  or  may  dissipate 
aloft  without  being  sensed  at  the  surface.  Either  of 
the  types  can  produce  dangerous  tailwind  and  downflow 
wind  shear  which  can  reduce  aircraft  performance  below 
critical  operating  levels,  especially  during  takeoff 
and  landing  (see  Fig.  1.1).  The  low-level  wind  shear 
associated  with  microburst  has  been  cited  as  a  main 
cause  of  many  a  -craft  accidents  in  the  U.S.;  for 
example,  see  studies  by  Fujita  and  Byers  (1977), 
Fujita  and  Caracena  (1977),  NT3S  (1983),  and  more 
recently  Fujita  (1985). 

The  above  definition  of  raicrobursts  was  somewhat 


Figure  l.l  .Schematic  illustration  of  an  aircraft 
penetration  of  a  aicroburst  (not  to  scale),  indic¬ 
ating  first  a  rapidly  increasing  headwind 
(perf ormance  increasing),  then  the  remnants  of  the 
downdraft  (performance  decreasing),  followed  by  a 
strong  tailwind  (seriously  performance 
decreasing),  AFTEP-  K* CAR-THy ^^7.  av 


the  aviation  involvement  (  McCarthy  e_t  a_l  .  ,  198  3). 
They  defined  a  microburst  as  a  downdraft-induced, 
damaging  horizontal  flow  near  the  surface,  whose  hori¬ 
zontal  dimension  is  <  4  ka,  and  whose  differential 
velocity  is  >  10  m  s*.  If  the  initial  horizontal 
dimension  is  >  4  km,  the  feature  is  referred  to  as  a 
macroburst  (Fujita  and  Wakimoto,  1983). 

In  an  effort  to  distinguish  microburst  types,  two 
subcategories  were  created.  According  to  Fujita  and 
Wakimoto  (1933),  a  "wet"  aicroburst  is  a  microburst 
event  also  accompanied  by  at  least  0.01"  of  precipita¬ 
tion  reaching  the  surface  between  the  onset  and  end  of 
the  high  winds,  including  any  calm  periods,  while  a 
"dry"  microburst  is  an  event  with  less  than  0.01"  of 
precipitation.  They  further  defined  the  term  "surface 
microburst"  as  a  microburst  with  its  outburst  winds 
detectable  by  a  network  of  ground-based  anemometers 
(Fig .  1.2). 

Mueller  and  Hildebrand  (1983)  detailed  the  struc¬ 
ture  of  a  raicroburst  and  compared  the  vertical  velo- 


(1983)  compared  the  low-level  horizontal  flow  computed 
from  dual-Doppler  radar  data  and  Portable  Automated 
Mesonet  (PAM)  surface  data.  Results  indicated  that 
the  derived  horizontal  velocities  were  faster  than 
those  observed  within  the  PAM  network.  They  attrib¬ 
uted  this  result  to  the  presence  of  frictional  forces 
in  the  surface  boundary-layer  which  were  at  the  level 
of  the  PAM  sensors  but  below  the  lowest  level  the 
radars  were  able  to  detect.  Wilson  _et_  a_l.  (  1984) 
identified  70  microbursts  using  Doppler  radar  results 
during  the  86  day  observation  period.  Using  PAM  data, 
Fujita  and  V a k i a o t o  (1983)  identified  up  to  136  micro- 
bursts  during  the  same  period.  At  this  point  it  is 
constructive  to  consider  what  is  known,  in  general, 
about  microbursts,  realizing  that  the  observed  samples 
from  which  this  climatology  is  drawn  is  not  only 
limited  in  size  but  in  geographical  extent  as  well. 

1.1  Microburst  Climatology 

During  the  last  decade,  two  extensive  field 
studies  have  been  conducted  observing  and  studying 
microbursts.  The  Northern  Illinois  Meteorological 
Research  On  Downbursts  (NIMROD)  was  conducted  in  the 
Chicago  area  from  19  May  through  1  July  1978,  lasting 
43  days.  The  Joint  Airport  Weather  Studies  Project 
(JAWS)  was  conducted  in  the  Denver  area  from  15  May 
through  9  August  1982,  totaling  86  days.  It  is  impor- 


tant  to  remember  when  comparing  results  of  these 
studies  that  they  occurred  in  areas  with  different 
environmental  conditions.  The  NIMROD  study  included 
the  effects  of  a  much  smoother  terrain  and  the  effects 
of  a  large  lake  (Lake  Michigan)  to  the  east,  while  the 
JAWS  environment  was  more  that  of  the  High  Plains  with 
little  locally  available  moisture  and  a  large  mountain 
range  immediately  to  the  west. 

Figures  1.3  and  1.4  show  the  number  of  micro¬ 
bursts  in  chronological  fashion  that  Fujita  and 
W'akimoto  (1983)  found  looking  at  the  PAM  data  from 
both  studies.  Notice  that  they  are  somewhat 
dissimilar.  The  NIMROD  data  appear  to  contain  many 
more  wet  occurrences  (32/50)  than  dry  ones  (13/50). 
On  the  other  hand,  the  JAWS  data  reveal  the  opposite. 
Among  the  186  microbursts  identified  by  PAM  data,  155 
were  dry  and  only  31  were  wet.  Moisture  availability 
and  the  higher  cloud  bases  generally  noted  in  the  High 
Plains  thunderstorm  allowing  precipitation  to  evapo¬ 
rate  prior  to  reaching  the  surface  could  account  for 
this  difference.  Also  note  the  apparent  cyclic  occur¬ 
rence  in  the  NIMROD  data.  Fujita  and  Wakimoto  (1983) 
suggested  that  this  may  be  due  to  synoptic  scale 
forcing  mechanisms  in  the  case  of  NIMROD  and  more 
local  effect  mechanisms,  such  as  orographic  lifting 
and  solar  heating  ,  in  the  case  of  JAWS. 


Figure  1.5  shows  the  occurrence  of  microbursts  by 
hour  for  the  two  studies.  The  diurnal  signal  seems  to 
be  present  in  the  JAWS  data  but  not  the  NIMROD  data. 
Also,  note  that  the  diurnal  signal  in  the  JAWS  data  is 
present  in  both  the  wet  and  dry  cases.  The  larger 
number  of  microbursts  which  occurred  after  daylight  in 
the  Midwest  can  be  attributed  to  the  larger  number  of 
nocturnal  thunderstorms.  Perhaps  the  most  revealing 
results  are  shown  in  Figs.  1.6  and  1.7.  Figure  1.6 
depicts  the  relationship  between  maximum  reflectivity 
in  the  parent  s'.orm  at  500  m  AGL  and  the  maximum 
differential  horizontal  wind  speed  across  the  micro¬ 
burst.  The  fact  that  the  correlation  is  almost  non¬ 
existent  suggests  to  a  forecaster  that  perhaps  he 
should  look  elsewhere  for  a  telltale  signal  of  an 
impending  microburst,  especially  when  making  a  wind 
shear  intensity  forecast  or  "nowcast".  The  rela¬ 
tionship  between  the  rainfall  rate  and  the  maximum 
wind  speed  is  shown  in  Fig.  1.7.  Again  the  correla¬ 
tion  is  very  small  suggesting  little  or  no  relation¬ 
ship  between  the  horizontal  wind  speed  and  the  amount 
of  rainfall  reaching  the  ground.  Figure  1.8  depicts 
front-rear  windshear  for  the  186  microbursts  identi¬ 
fied  by  Fujita  and  Wakimot.o  (1983).  While  one  may  be 
tempted  to  suggest  that  the  presence  of  precipitation 


enhances  the  horizontal  wind  shear,  the  reader  is 
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Ftg.  1.5  Diurnal  variation  of  surface  microburscs. 
Top:  186  aicroburscs  depicted  during  86  days  of 
JAWS  operation.  Two  peak  activity  periods  are 
1500  and  1800  MDT.  Bottom:  50  microburscs  during 
43  days  of  NIMROD  operation.  Rather  uniform 

activity  is  seen  from  1000  MDT  to  midnight. 
(_C.FT£E-  FUXi  T*\  ,  1963.  ) 
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Fig.  1.6  Correlogram  of  maximum  radar  reflec 
tivity  at  500  m  AGL  versus  maximum  velocity  dif 
ferential  (after  McCarthy  jst_  al_.  ,  1983). 
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Fig.  1.7  Scatter  diagram  of  the  relationship 
between  rainfall  race  and  the  maximum  wind  speed 
of  186  microbursts  during  JAWS 4c.) 


Fig.  1.6  Frequency  of  surface  aicrobursts  pre¬ 
sented  as  a  function  of  the  front-side  to  back¬ 
side  wind  shear.  An  aircraft  penetrating  through 
a  microburst  from  the  front  or  the  backside  would 
experience  this  magnitude  of  wind  shear. 


reminded  of  the  small  size  and  local  nature  of  this 
sample. 


Figure  1.9  shows  expected  changes  one  is  likely 
to  encounter  when  looking  at  macrobursts  and  micro- 
bursts.  For  a  microbursts  the  downflow  may  extend 
very  close  to  the  ground  prior  to  spreading  out.  In 
cases  where  the  precipitation  shaft  does  not  reach  the 
surface,  the  downward  moving  air  is  heated  dry-adiaba- 
tically  unless  evaporative  cooling  can  occur  (Fig. 
1.10).  When  precipitation  is  occurring  from  the  cloud 
base  to  the  surface,  cooling  is  less  likely  since  the 
air  is  already  highly  saturated.  Srivastava  (1935) 
notes  th3t  the  downward  moving  air  is  heated  nearly 
dry  adiabatically .  This  suggests  that  while  the  air 
is  descending  through  the  relatively  saturated  area 
within  the  precipitation  shaft  below  the  cloud  base, 
it  does  not  evaporate  enough  of  the  rain  to  either 
cool  the  air  or  cause  it  to  heat  up  at  the  moist 
adiabatic  rate,  possible  to  obtain  a  temperature 


Fig.  (-9  Three  stages  of  a  descending  microburst .  A  mid¬ 
air  microburst  may  or  may  not  descend  to  the  surface.  If 
it  does,  the  outburst  winds  develop  immediately  after  its 
touchdown .  (Arlene,  FwViTa  i<?3 S) 
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Fi g.  1-10  Schema tic  views  of  wet  and  dry  microbursts .  Wet  { 

microbursts  are  expected  to  occur  in  the  wet  regions  of  the  , 

world,  while  dry  m icrobursts  are  commonly  seen  in  the  dry  ' 

regions  with  high  bases  of  convective  clouds. (fibdtJZ  F=LxrrTA3/<)f5'J  i 


CHAPTER  2 


Statement  of  the  Problem 

The  main  objective  of  this  dissertation  research 
is  to  study  the  kinematic,  dynamic  and  thermodynamic 
(via  retrieval  techniques)  structures  of  a  microburst- 
producing  thunderstorm.  The  storm  which  produced  this 
nicroburst  was  of  "average"  strength.  One  of  the 
primary  reasons  a  storm  of  this  strength  was  chosen 
was  to  test  the  subtleness  of  the  techniques  in 
recovering  the  derived  fields.  Data  used  in  the 
research  are  dual-Doppler  radar  measurements  taken 
during  the  JAWS  project.  The  specific  storm  being 
investigated  occurred  on  14  July  1982  with  the  micro¬ 
burst  beginning  near  1641  MDT  and  continuing  through 
1651  MDT.  During  this  period,  a  peak  velocity 
difference  across  the  microburst  of  21  m  s*  was  ob¬ 
served,  see  McCarthy  e_t  al .  (  1983).  The  data  being 
considered  in  this  study  are  for  a  full  volume  scan 
beginning  at  1646  MDT  and  ending  at  1648  MDT. 

The  analysis  and  reduction  of  the  data  are  simi¬ 
lar  to  those  developed  by  Lin  and  Pasken  (1982)  and 
Wirsing  (1985).  Data  retrieved  will  include  horizontal 
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wind  fields  as  suggested  by  Armijo  (1969)  and  vertical 
wind  fields  using  the  anelastic  continuity  equation 
with  the  lower  boundary  condition  applied  similar  to 
Rodi  e_t  a_l..  (  1983).  Perturbation  pressure  and  temper¬ 
ature  deviations,  from  their  horizontal  averages,  will 
be  retrieved  using  the  method  outlined  by  Gal-Chen 
(1973).  These  data  will  be  judiciously  processed  and 
carefully  checked  using  the  objective  analysis  schemes 
developed  at  NCAR  and  here  at  St.  Louis  University. 
Once  the  thermodynamic  fields  are  retrieved,  momentum 
and  time  continuity  checks  proposed  by  Gal-Chen  and 
Hane  (1981)  will  be  employed  to  determine  the 
"goodness  "  of  the  derived  fields. 

Once  this  has  been  accompl ished ,  the  data  will  be 
used  to  test  the  hypothesis  that  a  downdraft  which 
produced  the  microburst  at  low  levels  was  triggered 
and  to  some  extent  maintained  by  a  favorable  perturba¬ 
tion  pressure  gradient  force  (PPGF)  together  with  the 
combined  effects  of  evaporative  cooling  and  precipita¬ 
tion  loading  in  the  middle  troposphere.  For  a  wet 
microburst,  precipitation  loading  is  more  important 
than  evaporative  cooling.  As  a  downdraft  approaches 
the  lower  troposphere,  a  strong  downflow  with  speeds 
reaching  -10  m  s1  impinges  upon  the  surface  forming 
a  "stagnation  mesohigh".  This  high  is  surrounded  by 
low  pressure  in  the  regions  with  the  strongest  outflow 
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resulting  in  a  pronounced  horizontal  PPGF  outward  from 
the  high  pressure  center.  This  PPGF  is  primarily 
responsible  for  driving  intense  diverging  outflow  in 
the  lower  layers,  extending  from  the  surface  to 
approximately  1  km  height. 

In  order  to  determine  the  role  played  by  friction 
in  thermodynamic  retrieval,  a  frictional  parameter¬ 
ization  scheme  similar  to  that  of  Klemp  and  Wilhelmson 
(1978)  will  be  employed.  The  variability  of  eddy 
viscosity  coefficients  and  Reynolds  stresses  (eddy 
transfer  of  momentum)  in  the  subcloud  layer  of  a  storm 
will  be  studied  in  detail.  Such  a  study  will  reveal 
the  effect  of  low-level  transfer  of  momentum  by  eddies 
on  the  structure  and  internal  dynamics  of  a  nonsevere 
thunderstorm  which  produced  a  microburst  in  Colorado. 


CHAPTER  3 

Data  Analysis  and  Reduction 

The  success  of  thermodynamic  retrieval,  to  a  con¬ 
siderable  extent,  depends  on  the  degree  of  accuracy  of 
the  detailed  wind  field  (Hane  e_t  al .  ,  1981;  Lin  and 

Pasken,  1982).  The  data  set  used  in  this  study  was 
obtained  from  the  JAWS  project.  It  represents  an 
individual  deep  volume  scan  of  a  storm  by  two  separate 
radars  beginning  at  1646  MDT .  The  storm  was  producing 
a  microburst  at  that  time  which  actually  lasted  from 
1643  to  1651  MDT. 

The  parent  storm  which  produced  the  microburst  is 
closer  to  an  "average"  strength  storm.  A  storm  of 
this  strength  is  more  likely  to  be  encountered  than 
storms  which  also  produce  severe  weather.  One  of  the 
primary  reasons  this  storm  was  chosen  for  study  was  to 
test  the  subtleness  of  the  techniques  in  recovering 
the  derived  fields  from  a  less  intense  storm.  In 
addition,  a  deep  volume  scan  is  available  for  this 
storm  from  the  surface  to  approximately  8  km  AGL. 
Hence  the  entire  storm  structure  can  be  investigated 
using  a  variational  approach  with  downward 
integration  . 

While  no  dynamic  studies  of  this  specific  storm 
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have  been  performed,  Wilson  e_t  aj^.  (1984)  studied  the 
kinematic  aspects  of  this  particular  microburst  at  the 
50  m  level  from  1641  to  1650  MDT.  In  addition, 
Mueller  and  Hildebrand  (1983)  analyzed  a  microburst 
case  from  the  JAWS  data  for  29  June  1982.  Results 
from  these  studies  will  be  compared  with  the  results 
from  the  present  study  in  Chapter  6. 

The  data  from  NCAR's  CP2  and  CP4  Doppler  radar 
were  chosen  for  use  in  this  study.  There  are  two 
reasons  for  choosing  these  radars.  First,  the  cross¬ 
beam  angle  formed  by  the  two  radars  is  larger  (78°  ) 
for  CP2  and  CP4  than  either  CP2  and  C?3  (43')  or  CP3 
and  C?4  (35°).  Second,  the  storm  and  the  radars 
nearly  form  an  isosceles  triangle  with  the  storm  at 
the  vertex  of  the  sides  of  equal  length.  These  tend 
to  minimize  the  effect  of  errors  and  variance  of  the 
measured  wind  field  on  the  derived  wind  field  (Fig. 
3.1). 

Studies  which  focus  on  the  kinematic  aspects  of 
storms  generally  require  only  first  derivatives  of 
derived  velocity.  This  study  utilizes  the  complete 
momentum  equations  and  their  differential  form  and, 
therefore,  requires  higher  order  derivatives.  Thus, 
the  criteria  used  in  the  processing  of  the  measured 
data  must  be  more  stringent  than  those  employed  in  a 
kinematic  study. 
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Section  3.1  describes  the  general  synoptic  situa¬ 
tion.  In  Section  3.2,  the  characteristics  of  the 
radars  and  specific  information  on  the  data  set  are 
detailed.  Section  3.3  discusses  the  data  reduction 
techniques  in  detail.  Each  critical  choice  is  re¬ 
viewed  and  reasons  for  specific  choices  offered.  The 
computational  procedure  to  synthesize  a  3-D  wind  field 
from  radial  velocities  and  radar  reflectivities  are 
presented  in  Section  3.4.  Section  3.5  discusses  the 
variational  adjustment  of  the  3-D  wind  field  to  con¬ 
form  to  the  lower  kinematic  boundary  condition  at  the 
surface.  Finally,  Section  3.6  discusses  the  recovery 
of  thermodynamic  vairables,  e.g.,  perturbation  pres¬ 
sure  deviation,  perturbation  density  deviation  and 
perturbation  virtual  tempera  ure  deviation,  in  detail. 

3.1  Synoptic  Situation 

The  synoptic  situation  at  0600  MDT  on  14  July 
1982  is  shown  in  Fig.  3.2.  It  revealed  a  high  pres¬ 
sure  center  in  northwestern  Colorado  with  a  cold  front 
farther  west.  This  cold  front  passed  through  the 
Denver  area  after  the  1800  MDT  14  July  time. 

The  sounding  for  Denver  released  at  0000  GMT  on 
15  July  (equivalent  to  1800  MDT  14  July)  is  presented 
in  Fig.  3.3.  It  shows  a  dry  layer  from  the  surface  to 
the  reported  cloud  base  at  3.5  km  AGL  (  Rodi  e_t_  al .  , 
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are  shown  in  Tables  3.1,  3.2,  and  3.3.  Different 
scanning  strategies  were  employed  due  to  joint  use  of 
CP2.  This  radar  was  required  to  make  a  360  degree 
sweep  at  low  elevation  angles  every  5  minutes.  This 
took  approximately  1  minute  to  complete  and  utilized 
the  capability  of  the  radar  to  change  pulse  repetition 
frequency  (PRF).  In  addition,  when  not  needed  by  the 
JAWS  personnel,  this  radar  made  a  complete  360  degree 
volume  scan  every  15  minutes.  The  other  two  radars, 
CP3  and  CP4,  were  used  for  JAWS  objectives. 

To  reduce  temporal  and  spatial  problems  within 
the  data  set,  an  effort  was  made  to  scan  the  same 
volume  at  the  same  time  with  more  than  one  radar.  The 
scanning  sequence,  including  the  elevation  angles 
selections  for  each  radar,  was  coordinated  between 
radars.  This  resulted  in  a  data  set  which  has  a  higher 
degree  of  temporal  coincidence  than  has  been 
previously  produced. 

Theoretically  and  ideally,  the  volume  scanned  by 
the  two  radars  should  be  the  same.  To  achieve  this 
with  two  Doppler  radars  can  be  somewhat  of  a  problem. 
It  involves  a  particular  geometric  interrelationship 
between  the  two  radars  and  the  storm.  The  storm 
should  be  located  so  as  to  form  an  isosceles  triangle 
with  the  storm  at  the  vertex  of  the  equal  legs. 
Figure  3.4  shows  the  relationship  between  the  radars 
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1933).  This  coincides  wich  the  decrease  in  temper¬ 
ature-dewpoint  spread  at  500  mb.  From  500  mb  the 
sounding  is  almost  moist  adiabatic  to  450  mb  where  a 
small  inversion  exists.  A  parcel  lifted  moist  adiaba- 
tically  from  530  mb  would  wipe  out  this  inversion  to 
280  mb.  Thus  with  very  little  forced  lifting,  a 
parcel  could  easily  rise  from  the  surface  to  270  mb 
(8.4  km  AGL).  Surface  temperature  was  30  C  and  the 
lapse  rate  was  nearly  dry  adiabatic  from  the  surface 
to  600  mb.  This  sounding  is  similar  to  the  composite 
microburst  sounding  presented  by  Brown  _e_t  a_l  .  (  193  2). 

The  Program  for  Regional  Observing  and  Forecas¬ 
ting  Service  (PROFS)  noted  early  in  the  day  the  poten¬ 
tial  for  the  production  of  dry-type  microbursts.  This 
proved  to  be  the  case  with  7  microbursts  identified 
during  that  day.  The  dry  mid-level  and  low-level 
conditions  were  conducive  to  the  enhancement  of  down¬ 
ward  moving  air  when  evaporative  cooling  in  the  sub¬ 
cloud  region  increased  negative  buoyancy  (Srivastava, 
1985)  . 

3.2  Doppler  Measurements 

The  data  used  in  this  study  were  obtained  from  two 
Doppler  radars  used  in  the  JAWS  project.  These 
radars,  designated  CP2  and  CP4,  were  located  in  the 
area  of  Stapleton  International  Airport  in  Denver,  CO. 
(Fig.  3.4).  The  technical  parameters  of  these  radars 


and  the  microburst  of  this  study.  CP2  and  CP4  nearly 
form  the  ideal  configuration  to  take  advantage  of  this 
situation.  The  more  orthogonal  the  radar-storm-radar 
angle  the  less  the  effect  of  variance  of  the  sensed 
radial  wind  on  the  possible  range  of  the  true  wind. 
Figure  3.1  shows  that  a  wider  range  of  true  winds  is 
possible  when  the  angle  between  the  storm  and  the 
radars  is  decreased. 

In  both  cases  (  Fig.  3.1)  the  two  radars  sense 
the  same  vector  wind  plus  or  minus  the  same  variance. 
Any  vector  wind  originating  at  point  P  and  terminating 
within  the  outlined  area  will  produce  the  correct 
sensed  radial  wind  plus  or  minus  the  variance.  It  is 
evident  from  the  size  and  shape  of  the  areas  that  when 
the  angle  between  the  radars  and  the  storm  is  90 
degrees  the  size  of  the  area  is  minimal. 

3.3  Data  Analysis  Procedures 

The  Doppler  data  were  brought  back  from  NCAR 
having  been  edited  as  described  in  the  previous  sec¬ 
tion.  They  consisted  of  radial  velocity,  standard  de¬ 
viations  of  radial  velocity  (referred  to  as  spectral 
width  hereafter),  and  reflectivities.  The  spectral 
widths  were  of  no  immediate  use  due  to  the  random 
errors  introduced  during  the  collection  procedure 
(Elmore,  1984;  private  communication). 


Lin  and  Pasken  (1982)  discussed  six  stages  in 
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analysis  of  the  Doppler  data.  They  are: 

1)  Desired  data  are  collected  and  ordered  accord 
ing  to  time  and  altitude. 

2)  Errors  are  deleted  and  aliased  velocities,  or 
"folded  velocities",  are  corrected. 

3)  Data  are  then  interpolated  on  to  the  three- 
dimensional  grid  array. 

4)  The  three-dimensional  wind  field  is  computed 
from  dual-Doppler  radar  data;  u,  v,  and  w  components 
are  then  calculated  from  gridded  fields  of  reflec¬ 
tivity  and  radial  velocity. 

5)  Kinematic  variables  are  derived,  and 

6)  These  data  are  displayed. 

The  specific  detail  of  each  of  these  steps  will 
be  covered  with  the  rationel  for  the  choices  of  tech¬ 
nique  used  in  this  study. 

3.3.1  Extraction  of  desired  parameters 

The  data  brought  back  to  St.  Louis  University 
consisted  of  a  single  magnetic  tape  for  each  radar. 
The  first  step  was  to  reduce  this  vast  quantity  of 
data  into  a  more  compact  and  usable  form  for  the  PDP 
11/70  computer.  The  data  were  first  scanned  and  any 
data  falling  outside  our  predefined  area  of  interest 


were  eliminated. 


Next  the  data  were  converted  from 
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spherical  coordinates  centered  on  each  radar  to  an 
orthogonal  coordinate  system  centered  on  radar  CP2. 

During  this  step  the  location  of  the  data  was  trans¬ 
lated  using  storm  motion  and  actual  observation  time 
of  the  individual  ray  scans.  Finally,  the  data  were 
"cut"  horizontally  in  relationship  to  the  desired 
analysis  levels.  A  slab,  0.5  km  thick,  centered  on 
each  of  the  analysis  levels  was  defined  and  all  obser¬ 
vations  within  that  slab  were  consolidated  into  a 
single  file.  The  choice  of  .5  km  was  dictated  by  need 
to  have  enough  data  points  to  satisfy  the  Barnes 
criteria  for  minimum  number  of  observations  (Barnes, 

1973).  This  file  was  further  reduced  by  limiting  the  j 

actual  areal  coverage  of  the  observations  (Fig.  3.5).  1 

i 

i 
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3.3.2  Quality  control  of  extracted  parameters  j 

It  is  important  that  the  data  be  rendered  con-  1 

sistent  with  the  known  features  of  the  storm.  The  j 

quality  control  checks  used  to  insure  this  consistency 
were  as  follows: 

1)  All  reflectivities  were  checked  for  values  [ 

greater  than  20  dBz.  This  value  corresponds  well  to 

the  cloud  boundary.  In  cases  where  the  reflectivity 
was  less,  the  observation  was  flagged  as  suspect,  but 
not  eliminated  from  the  data  base. 

2)  All  velocities  were  again  checked  for  folding. 

This  involved  objectively  computing  radial  and  cross- 


beam  shear  and  comparing  these  to  maximum  allowable 
values.  In  cases  where  the  shear  was  excessive,  the 


suspect  data  were  unfolded.  Finally  the  unfolded 
velocities  were  subjectively  checked  in  areas  which 
appeared  suspicious. 

3)  At  this  point  no  quality  control  checks  were 
considered  for  the  spectral  widths. 

4)  Finally,  the  velocities  were  checked  in 
conjuction  with  reflectivities  and  the  apriori  know¬ 
ledge  of  storms  of  this  type  and  the  environmental 
conditions  to  insure  against  gross  errors  which  might 
be  caused  by  non-hydrometer  targets  within  the  radar 
scan,  such  as  birds,  insects,  planes,  etc. 

3.3.3  Interpolation  of  the  extracted  parameters  to 
grid 

The  actual  steps  required  to  analyze  the  data 
will  be  presented  here.  In  order  to  obtain  as  consis¬ 
tent  a  set  of  data  as  possible,  at  each  step  the  data 
were  subjected  to  both  objective  error  analysis  and 
scrutiny . 


3. 3. 3.1  Collection  and  ordering  of  desired  data 

This  procedure  was  basically  completed  with  the 
cutting  of  the  data  in  the  horizontal.  However,  a 
correction  was  applied  to  account  for  the  time  differ¬ 
ence  between  the  different  sweeps  which  made  up  the 


horizontal  slab  of  observations.  The  entire  volume 
scan  for  both  radars  began  at  16:46:15  MDT  and  ended 
at  16:43:15  MDT. 

Wilson  (1933)  suggested  a  temporal  correction 
scheme  which  also  accounted  for  the  fact  that  the 
radial  velocity  is  a  vector  quantity  and  not  a  scalar. 
As  such,  any  correction  in  space  would  also  require  a 
realignment  of  the  vector  to  maintain  its  length  and 
sense.  In  situations  where  the  elevation  angles  of 
the  radar  are  high,  this  can  result  in  an  error  in  the 
corrected  velocities.  However,  in  this  study,  the 
maximum  horizontal  wind  error  would  be  less  than  1.0 
m  s  . 

Parsons  et_  a_l_.  (  1933)  employed  an  advection 
correction  scheme  suggested  by  Gal-Chen  (1932)  to 
account  for  the  temporal  differences  in  scanning  the 
storm  volume.  Vertical  velocity  errors  were  reduced 
by  an  order  of  magnitude.  They  noted  that  the  error 
on  computed  vertical  velocities  was  a  function  of 
three  parameters;  namely,  the  scanning  time,  the 
advection  velocity,  and  the  correction  velocity. 
Increased  scanning  time  increased  the  errors  in  verti¬ 
cal  velocity  both  in  the  corrected  and  uncorrected 
cases.  It  also  threatens  the  assumption  of  steady 
state  used  in  studies  involving  only  one  time  scan. 
Similar  results  also  held  for  the  advection  velocity. 


The  errors  in  vertical  velocity  seemed  to  be  ne g- 
atively  correlated  with  corrected  speed,  i.e.,  higher 
vertical  velocities  produced  more  negative  errors. 
However,  errors  were  still  much  smaller  for  corrected 
cases  when  compared  to  uncorrected  cases.  They 
concluded  that  "...  considerable  improvement  can  be 
obtained  in  the  derived  wind  fields  even  when  the 
advection  velocity  is  not  exactly  known." 

The  advection  scheme  was  employed  to  correct  the 
data  location.  The  difference  between  the  "storm 
observation  time"  and  the  actual  data  observation  was 
multiplied  by  the  mean  storm  motion  to  produce  a 
horizontal  offset.  This  was  then  added  to  the  hor¬ 
izontal  coordinate  to  produce  a  new  location  for  the 
data  point,  so  that 

Xq  =  x0-  ( t  -  t0  )  cx  (3.3.1) 

yq=y<k-(t-to)c3r  (3.3.2) 

where 

( x0  ,  y0 )  a  uncorrected  datura  coordinates, 

( xq  ,  yq  )  *  corrected  datum  coordinates, 

(cx.cy)  =  mean  storm  velocity  in  x  and  y  direction, 
tD  a  mean  storm  observation  time  (averaged  over 
the  starting  and  ending  times),  and 
t  a  actual  radar  time  for  datura  point. 

3. 3. 3. 2  Interpolation  scheme 
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There  are  a  number  of  techniques  available  for 
transferring  the  data  to  a  grid.  Lhermitte  (1970) 
utilized  a  COPLAH  technique  in  a  cylindrical 
coordianate  system.  The  advantage  of  this  system  was 
that  Cartesian  wind  components  could  be  synthesized 
directly  from  cylindrical  components.  By  contrast, 
many  of  the  studies  involving  JAWS  data  use  the  Car¬ 
tesian  grid,  e.g.,  Mueller  and  Hildebrand  (1983), 
Kessinger  e_t_  a_l. •  (1983),  Rodi  e_t_  al_ .  (1983),  etc.  In 
keeping  with  this  precedent,  the  Cartesian  grid  system 
is  utilized  in  this  study.  There  are  a  number  of 
drawbacks  to  this  choice.  The  first  involves  the 
necessity  to  convert  radial  velocities  from  two  radars 
to  horizontal  velocities.  This  involves  the  solution 
of  a  hyperbolic  partial  differential  equation  to  de¬ 
rive  the  vertical  wind  (Armijo,  1969).  The  second  and 
perhaps  more  bothersome  problem  involves  transferring 
winds  collected  on  a  conical  surface  to  a  horizontal 
plane  . 

When  the  data  are  cut  into  slabs  the  location  of 
data  points  within  that  slab  is  not  random.  As  a 
result,  the  technique  for  evaluating  the  variable  at 
grid  points  must  be  carefully  considered.  There  are 
many  techniques  for  evaluating  non-gridded  data  at 
grid  points.  Only  those  which  involve  some  averaging 


process  of  the  points  "close"  to  the  grid  point  will 
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be  considered. 

The  simplest  of  these  involves  a  direct  unweighted 
average  of  all  points  within  a  specified  radius.  This 
technique  is  inappropriate  since  it  smooths  the  data 
far  too  much,  removing  many  of  the  mesoscale  features 
of  interest.  Cressman  (1959)  suggested  a  scheme  which 
weighted  the  influence  of  the  data  in  relation  to  the 
distance  from  the  grid  point: 


where  R  is  the  maximum  radius  of  influence  and  r  is 
the  actual  distance  from  the  grid  point  to  the  datum 
point.  This  scheme  has  been  used  by  Mueller  and 
Hildebrand  (1983),  Kessinger  e_t  al_.  (1983)  and  Wilson 
et  al .  (1984).  The  primary  criticism  of  this  techni¬ 
que  is  the  large  amount  of  influence  this  technique 
exerts  on  a  grid  point  at  great  distances  and  the 
requirement  from  multiple  passes  at  different  radii  of 
influence  to  reintroduce  the  effects  of  smaller  scale 
features.  Perhaps  even  a  more  pointed  criticism  is 
the  inability  to  produce  a  "response  curve"  similar  to 
that  found  with  the  Barnes  technique.  By  choosing 
differing  radii,  the  actual  amplitude  recovered  from  a 
specific  wave  number  is  not  readily  apparent.  Thus, 
the  investigator  must  be  careful  when  discussing  fea¬ 
tures  which  may  be  below  the  range  of  discernible  wave 
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numbers . 

Barnes  (1964)  developed  a  weighted  influence 
function  which  partially  negated  the  problems  of  the 
Cressman  scheme.  The  result  was  an  exponential 
weighting  function  which  could  be  reapplied  succes¬ 
sively  to  obtain  a  gridded  set  of  data.  It  still, 
however,  required  the  use  of  multiple  passes,  i.e., 


where 


1-  e*r 


(3.3.4) 


R  is  the  radius  of  influence, 

r  is  the  distance  between  the  datum  point  and  the 
grid  point ,  and 

E  is  the  amount  of  influence  of  the  weighted  data 
we  wish  to  represent. 

When  E  =  4.0,  it  implies  that  we  have  represented  98 
percent  of  the  influence  of  any  datum  within  a  circu¬ 
lar  region  whose  radius  is  R.  This  scheme  was  im¬ 
proved  by  Barnes  (1973)  reducing  the  number  of  passes 
to  2,  namely. 


e*p 


(3.3.5) 


where  R,  r,  and  E  are  as  above  and  is  a  parameter 
between  0  and  1  which  represents  confidence  in  the 
data  (Koch  £t_  a_l .  ,  1984). 

As  stated  above,  one  of  the  advantages  of  the 
Barnes  scheme  is  the  ability  to  compute  the  scheme's 
response  to  various  wave  number  inclusions  in  the 
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final  gridded  data  set.  Equation  (3.3.6) 
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(3.3.6) 


shows  the  computational  form  of  the  response  curve  an 
Fig.  3.6  shows  the  final  response  curve  for  the  two 
scan  radii,  R  =  1.25  km  and  R  =  1.75  km,  used  in  this 
study . 

In  most  of  the  studies  in  the  literature,  the 
Barnes  scheme  has  been  applied  primarily  to  the  hori¬ 
zontal  plane  with  either  the  p  or  z  coordinate  held 
constant,  e.g.,  Lin  and  Pasken  (19S2),  Wirsing  (1935), 
and  Lin  ejt  aj^ .  (1935).  An  in-depth  discussion  of  the 
Barnes'  gridding  technique  as  it  was  used  in  this 
study  can  be  found  in  Appendix  A. 

Employment  of  the  horizontal  gridding  scheme  re¬ 
quires  that  all  data  points  within  the  slab  be 
projected  directly  to  the  level  of  interest.  Since 
the  data  are  not  randomly  located  within  the  slab  but 
rather  in  intersection  bands,  the  resultant  wind  field 
contains  features  which  appear  as  wave  patterns. 
These  patterns  are  fictitious  and  are  the  result  of 
speed  and  directional  wind  shear  in  the  vertical. 
Figure  3.7  shows  this  phenomenon.  Point  one  obtains 
all  its  influence  from  the  points  to  the  right  which 
in  actuality  lie  above  the  plane,  while  the  point  on 


the  left  obtains  all  its  influence  from  the  points  to 
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PESPGNSE  CURVE  FOR  BARNES  SCHEME 


Fig.  3.6  Final  response  curve  from  Barnes' 
gridding  scheme  for  R  *  1.25  km  and  1.75  km. 
Gamma  ■  0.3.  Ordinate  represents  the  wavelength 
and  abscissa  is  the  percentage  of  original 
amplitude  retained  after  the  gridding  process. 
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Che  left  which  lie  below  the  plane.  The  result  is 
that  the  vertical  directional  and  speed  shear  is 
transferred  to  the  horizontal  plane. 


Extension  of  the  Barnes  scheme  to  three  dimen¬ 
sions  has  not  appeared  extensively  in  the  literature 
and  where  it  has  there  is  considerable  variation  in 
the  method  for  handling  the  vertical  coordinate. 
Regan  (1984;  personal  communication)  utilized  a 
weighting  function  which  employs  a  spherical  volume  of 
influence.  Wang  (1934;  personal  communication) 
utilized  an  elliptical  volume  in  which  the  vertical 
axis  is  smaller  than  the  horizontal  axis  and  the 
horizontal  axes  are  equal.  Barber  (1983)  suggests  a 
weighting  function  of  the  form 


r 

( 


\ 


A* I  >.'s 


ev 


(3.3.7) 


where 

k  and  D  are  weight  parameters  controlling  hori¬ 
zontal  and  vertical  smoothness,  respectively, 

and 

d‘L  and  Z  ^  are  horizontal  and  vertical  distances, 
respectively,  between  observation  and  a  given 
gridpoint . 

This  essentially  results  in  an  elliptical  weighted 
volume,  with  K  and  D  defining  the  major  (horizontal) 


and  minor  (vertical)  axis,  respectively.  On  the  other 
hand,  Wirsing  (1985)  attempted  to  solve  the  problem  in 


a  completely  different  manner.  Along  with  gridding 
the  data  in  a  manner  similar  to  the  two-dimensional 
scheme,  he  also  gridded  the  heights  of  the  actual  data 
points.  The  gridded  values  were  then  assigned  to  the 
gridded  height  and  the  values  at  the  level  of  interest 
were  interpolated  in  the  vertical  from  the  gridded 
profiles.  This  technique  eliminates  the  requirement 
to  determine  whether  a  datum  point  is  within  an  ellip¬ 
tical  volume  and  adds  only  one  additional  run  through 
the  gridding  routine  and  one  additional  storage  array. 

This  technique  was  employed  in  the  present  study.  A 
number  of  different  fitting  polynomials  were  studied 
and  the  linear  fit  seems  to  be  the  most  sensible.  ’v'e 
feel  that  the  order  of  the  polynomial  used  should  be 
•'nversely  proportional  to  the  standard  deviation  of 
the  interval  of  adjacent  data  points  for  the  inde¬ 
pendent  variable.  Since  the  intervals  are  of  non- 

« 

| 

uniform  size,  the  simplest  polynomial  provides  the  j 

most  reliable  approximation. 

3.4  Synthesis  of  three-dimensional  wind  field 


The  method  of  processing  Doppler  data  to  obtain  3- 
dimensional  wind  fields  and  thermodynamic  information 
are  separate  procedures  in  this  study,  but  the  re- 


covery  procedures  for  the  thermodynamic  variables 
depend  on  the  procedures  used  to  derive  the  wind 
field.  Lhermitte  (1963)  and  Armijo  (1960)  first  sug¬ 
gested  methods  of  recovering  the  3-dimensional  wind 
field  patterns  from  multiple  Doppler  radar  data. 
Lhermitte  (1968)  developed  a  technique  referred  to  as 
the  COPLAN  scheme.  This  scheme  utilizes  cylindrical 
coordinates  and  processes  the  wind  field  on  tilted 
planes.  The  principle  drawback  of  the  COPLAN  scheme 
seems  to  be  the  requirement  that  the  volume  be  scanned 
in  unison  by  two  radars.  Zrnic  and  Doviak  (1934) 
provided  a  condensed  discussion  of  this  technique. 
They  noted  that  it  is  easier  to  compute  the  vertical 
velocity  using  the  COPLAN  technique  than  any  other 
principle  method  available,  e.g.,  Armijo  (1969). 

The  Armijo  (1969)  technique  for  2  Doppler  radars 
involves  direct  synthesis  of  the  3-dimensional  wind 
field  in  Cartesian  coordinates.  The  principle  draw¬ 
back  is  the  need  to  solve  an  inhomogeneous  hyperbolic 
partial  differential  equation  to  obtain  vertical  wind. 
However,  it  does  not  have  the  strict  synchronous  scan 
requirement,  rather  relying  on  a  temporal  adjustment 
scheme  to  adjust  for  nonsynchronous  scanning.  Since 
most  of  the  analysis  presented  in  this  study  employed 
the  interpolation  scheme  in  a  Cartesian  coordinate 
system,  the  Armijo  (1969)  method  of  producing  the  3- 


dimensional  wind  field  is  used. 


Considering  the  geometry  of  two  radars  looking  at 
the  same  storm,  Armijo  (1969)  expressed  the  radial 
velocities  as  sensed  by  the  radars  as  a  function  of 
the  three-dimensional  wind  field  by 


\r.,=  R,  (3.4.1a) 

v/  u(x-xj  -v(v-Vzi  -(W4Vj(iE-Zz) 

V^=  - — -  (3.4.1b) 

where  z- 

Rjis  racial  distance  from  radar  1  to  target, 

32i-s  radial  distance  from  radar  2  to  target, 

|  ^  s  r  3  d  i  a  1  velocity  measured  by  radar  1, 
is  radial  velocity  measured  by  radar  2, 

(u,v,w)  are  components  of  true  wind  in  x,y,z 

directions,  respectively, 

(x,y,z)  is  actual  location  of  the  target  w.r.t. 
origin , 

(xj  .  y,  i  z,  )  is  location  of  radar  1  w.r.t.  origin, 

( x2.'  *  z2, )  is  location  of  radar  2  w.r.t.  origin, 

and 

mean  terminal  velocity  of  hydrometeors 

within  the  target  volume. 

With  the  location  of  the  radars  and  the  measured 


velocities  known,  (3.4.1)  results  in  2 
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equations  in  4  unknowns,  namely,  u,  v,  v,  and  V^_  .  In 
order  to  solve  for  u,  v,  and  w,  it  is  necessary  to 
formulate  two  r-’ditional  equations.  Using  the  reflec¬ 
tivity,  a  terminal  velocity-reflectivity  equation  can 
be  formulated  relating  the  terminal  velocity  of  the 
target,  ,  to  reflectivity,  Ze.  A  number  of  inves¬ 
tigators  have  developed  such  emperical  formulae  in¬ 
cluding  Rogers  (1964),  Sekhon  and  Srivastava  (1971), 
Atlas  e_t  al .  (1973),  Battan  (1973),  and  Martner 
(1975).  Foote  and  du  Toit  (1969)  proposed  a  density 
height  correction  to  account  for  differences  in 
terminal  fall  speed  of  specific  sized  particles  due  to 
differences  in  air  density  through  which  the  particles 
are  falling.  The  Martner  formula  which  was  derived 
for  High  Plains  thunderstorms  along  the  Foote  and  du 
Toit  correction  are  employed  in  this  study,  i.e., 


where 


(3.4.2) 


is  raindrop  terminal  velocity, 

Ze  is  equivalent  radar  reflectivity  factor, 

is  standard  surface  density,  and 
o 

p  is  actual  environmental  density. 

This  formula  has  appeal  beyond  the  fact  that  it 
was  developed  for  High  Plains  thunderstorms.  In 


studies  where  the  reflectivity  is  small  or  may  be 


questionable,  a  form  of  this  equation  which  is  rel¬ 
atively  insensitive  to  errors  in  Ze  is  desirable. 
Formulas  which  utilize  small  exponents  tend  to  be  less 
sensitive  than  those  which  utilize  large  exponents. 
Table  3.4  compares  the  values  for  a  number  of 
equations  and  a  variety  of  reflectivities.  The  per¬ 
centage  change  from  10-60  dBz  reveals  that  while  the 
values  are  absolutely  larger,  the  percentage  change  is 
smaller.  Thus  the  Martner's  formula  is  chosen  for 
this  study. 

The  final  equation  is  the  anelastic  continuity 
equation : 


where 


02. 


(3.4.3) 


C(£) 


VI 


is  adiabatic-hydrostatic  density,  a  function 


of  height  only,  and 
is  horizontal  velocity  vector. 

Equations  (3.4.1),  (3.4.2),  and  (3.4.3)  form  a 
complete  set  of  equations  with  four  unknowns,  u,  v,  w 


and  V, 


The  solution  can  be  obtained  from  a  standard 


iteration  procedure  from  one  level  to  another.  The 
vertical  wind,  w,  is  consistent  with  the  horizontal 
winds  through  the  use  of  the  anelastic  continuity 
equation.  In  practice,  these  equations  are  modified 
and  solved  for  u,  v,  and  w.  Section  3.4.1  outlines 


this  procedure  in  detail. 


The  anelastic  continuity  equation,  Eq.  (3.4.3), 
requires  the  adiabatic  density  9<a.  •  Beginning  with 
the  first  law  of  thermodynamics 

dQ  =  Cp  dT  -  SdP  (3.4.4) 

where 

ft  is  specific  volume, 

Cp  is  specific  heat  at  constant  pressure, 
dT  is  change  in  temperature, 
dP  is  change  in  pressure, 

dQ  is  amount  of  internal  heat  added  to  or  lost 
from  the  system,  and  the  overbar  indicates 

environmental  (basic)  state. 

From  the  adiabatic  assumption  dQ  =  0,  one  has 

0  =  cpdT  -  adP  (3.4.5) 

The  hydrostatic  equation  gives 

dP  =■  -fg  (3.4.6) 

dz 

Substituting  Eq .  (3.4.6)  into  Eq .  (3.4.5)  yields 

dT  -  -  (3.4.7) 

Tz  Cp 

Replacing  temperature  using  the  equation  of  state 
1  dP  *  2  dP  -  g 

dz  j  dz  Tf  (3.4.8) 

Rearranging  Eq .  (3.4.8)  yields 

-I  ft  "Z 


V 


(3.4.9) 
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Integrating  this  equation  from  to  Z^.+-|  and  noting 
the  value  of  g  *  9.8  m/s,  we  arrive  at 

/- ZA.3ZO? 

p.  =  Pexp  — : -  (3.4.10) 

itvi  -u  1  \  T  u. 


t 


where  T^l+j  is  the  mean  environmental  virtual  temper¬ 
ature  between  layers  i  and  i+1.  To  solve  for  ^  >  we 
first  compute  the  value  of  from  the  surface  obser¬ 
vation  taken  at  1800  MDT ,  14  July  1982  in  Denver,  CO. 

To  determine  the  validity  of  this  technique  for 
computing  density  with  the  assumptions  made,  values  of 
density  were  computed  using  Eq .  (3.4.10)  and  the  equa¬ 
tion  of  state.  Table  3.5  shows  these  computations  as 
well  as  the  correction  factor  of  Foote  and  du  Toit 
using  the  different  density  values. 


3.4.1  Synthesis  of  wind  components 

To  solve  for  the  wind  components  we  begin  with 
Eq.  (3.4.1).  Multiplying  (3.4.1a)  by  (x-xz  )  and 
(3.4.1b)  by  (x-x,)  and  rearranging,  we  have 


Similarly  multiplying  (3.4.1a)  by  (y-y^)  aad  (3.4.1b) 
by  (y-y, ),  we  can  obtain 


-v-(v.v*V^ 


*V.  -  a*,  ±yz 


(3.4.11b) 


Table  3.5 


Comparison  of  Density  Computation  Methods 


Pres 


F  &  D 


F  &  D 


mb 

kg/m5 

CO 

835 

0.954 

1.096 

1.013 

1 .067 

700 

O.S43 

1.151 

0.936 

1 .104 

500 

0.658 

1.271 

0.853 

1.145 

400 

0.553 

1.363 

0.775 

1 .191 

300 

0.440 

1.493 

0.699 

1  .  240 

250 

0.384 

1.576 

0 . 628 

1.295 

0.320 

1.695 

0.562 

1.354 

60 


where 


=  y  - 


and  i 


1  or  2 


Equations  (3.4.11a)  and  (3.4.11b)  replace  those  of 
(3.4.1)  and  along  with  (3.4.3)  and  (3.4.4)  form  a  com¬ 
plete  set  of  equations. 

To  obtain  the  vertical  velocity  by  downward  inte¬ 
gration,  the  kinematic  boundary  condition  requires 
that  w  =  0  at  the  top  boundary.  This  height  was 
chosen  by  assuming  that  the  storm  top  could  be  defined 
by  the  20  d3Z  envelope.  At  the  level  actually  chosen 
as  the  storm  top,  Z^=S . 5  km,  the  maximum  reflectivity 
was  less  than  20  d3Z.  To  avoid  the  concern  that  per¬ 
haps  the  vertical  velocity  at  the  storm  top  was  not  in 
fact  zero,  a  cap  was  constructed  at  0.5  km  above  the 
highest  point  in  the  storm.  It  was  assumed  that  w  was 
actually  zero  at  this  height,  9.0  km.  Next  the  diver¬ 
gence  for  the  storm  top  level  is  comp"ted  from 
(3.4.11a)  and  (3.4. lib)  assuming  that  w=0  for  this 
level.  A  solution  for  vertical  velocity  is  obtained 
for  the  second  level  below  the  storm  top,  at  3.0  km 
Using  these  values  of  vertical  velocity  and  the  fact 
that  the  vertical  velocity  is  actually  zero  at  9.0  km, 
the  "second  guess"  of  the  vertical  velocity  is  made 
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for  the  storm  top  level  by  averaging  the  vertical 
velocities  from  9.0  and  3.0  km.  This  is  then  compared 
to  the  "first  guess"  of  the  velocities  at  8.5  km.  This 
iterative  process  is  continued  until  the  "nth  guess"  - 
"n-lth  guess"  is  within  a  specified  tolerance  level 
for  all  grid  points.  Rodi  e_t  a_l.  (  1983)  used  this 
scheme  with  upward  integration  technique  iterating 
until  the  variation  in  divergence  is  below  a  threshold 
value.  However,  they  continue  to  utilize  this  pro¬ 
cedure  throughout  the  volume.  Lin  and  Pas ken 
(1932a, b),  and  Wirsing  (1935)  applied  a  similar 
technique  except  they  began  at  the  top  of  the  storm 
and  integrated  downward.  The  above  technique  was  only 
used  to  obtain  the  vertical  velocity  at  the  topmost 
level.  From  that  point,  standard  downward  integration 
procedures  were  used. 

At  this  point  it  is  germane  to  address  the  use  of 
upward  integration.  Many  investigators,  e.g.,  Lin  and 
Pasken  (19S2a,b),  Ray  _e_t_  a_l_.  (  1980),  Nelson  (1980)  and 
Brown  and  Nelson  (1982),  suggested  that  downward  inte¬ 
gration  from  the  top  of  a  storm  volume  will  produce 
errors  of  less  magnitude  than  upward  integration  from 
the  bottom  of  the  storm  volume. 

3y  integrating  the  anelastic  continuity  equation 
to  obtain  the  vertical  velocity  at  the  next  level,  one 
can  see  that  errors  computed  early  in  the  downward 


integration  process  are  initially  multiplied  by 


smaller  density  terms.  While  integrating  upward, 
these  errors  would  be  influenced  by  the  larger  density 
terms  near  the  surface,  i.e., 

pu) , 

4r^-wi  ■*-  (3.4.12) 

where 


is  vertical  velocity  at  upper  level, 

W.  is  vertical  velocity  at  lower  level, 

is  horizontal  divergence  averaged  vertically 
through  the  layer  of  interest, 

■P(z)  is  basic  density  at  lower  level, 

JQj 

pa  U)  is  basic  density  at  upper  level, 

^2((z)  is  mean  density  averaged  vertically  between  two 
levels . 

In  this  study,  data  were  available  through  the 
entire  volume  of  the  storm  and  downward  integration 
was  used. 


3.5  Correction  of  Wind  Components  using  Variational 
Analysis 

As  indicated  previously,  computation  of  vertical 
velocity  was  accomplished  using  downward  integration 
of  the  anelastic  continuity  equation.  This  technique 
minimizes  errors  due  to  various  computational 


considerations  such  as  derivative  formulation  and 
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boundary  considerations,  and  errors  in  data.  Errors 
in  the  derived  vertical  velocity  are  generally 
greatest  near  the  surface  since  they  tend  to  be 
cumulative.  Techniques  which  correct  these  errors  are 
generally  based  on  a  linear  distribution  of  the 
surface  error.  The  method  suggested  by  O'Brien  (1970) 
is  one  such  technique  available  in  which  the  error 
variance  is  assumed  to  be  a  linear  function  of  height 
(or  pressure).  Another  scheme  suggested  by  Sasaki 
(1970)  involves  the  use  of  the  calculus  of  variation. 

Ziegler  e_t  a_l  .  (19S3)  employed  this  technique  in 
a  study  using  two  and  multiple  radars.  A  more  com¬ 
plete  discussion  of  this  technique  and  the  formulation 
of  the  Euler-Lagrange  equations  is  given  in  Appendix 
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weighting  factors  in  the  following  manner: 

C/rUerf'  •  (3.5.2) 

The  actual  computation  of  (J^  and  were 

similar  to  computation  of  the  horizontal  wind  com¬ 
ponent  equations  of  Armijo  (1969)  with  (C  set  to  zero. 
Wilson  _e_t  a_l.  (1984)  suggested  another  method  which 
seems  to  completely  remove  any  dependence  on  the 
sensed  data  and  is  related  only  to  the  relative  lo¬ 
cation  of  the  radar  and  an  initializing  assumption 

r'1  r"'"2" 

about  the  values  of  and  0^ .  This  technique  is 
employed  in  this  study  and  is  further  discussed  in 
Appendix  C. 

The  equations  used  in  the  adjustment  scheme  were 


-  -r  rT 

U=  U  -  f  UU 

v  =  V  -Hf 'wy  -Ju" 

(  V  *(  fv3)  -  0 


(3.5.3) 

(3.5.4) 

(3.5.5) 


Equations  (3.5.3)  and  (3.5.4)  form  the  Euler-Lagrange 
equations  used  in  the  variational  correction  tech¬ 
nique.  Equation  (3.5.5)  is  the  anelastic  continuity 
equation  which  acts  to  constrain  the  solution.  Using 
equations  (3.5.3)  and  (3.5.4),  we  arrive  at  the 


integral  elliptical  equation 

,0 


>jri?=  (f  ( D-S> ) at 


(3.5.6) 


■j  ,~w.  r --  j.-j.  • 
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where  D  is  the  "true"  horizontal  divergence 

/N 

D  is  the  computed  horizontal 
divergence 

A. 

and  D-D  is  the  horizontal  divergence  error 

Equation  (3.5.6)  represents  a  Poisson  type  partial 
differential  equation  which  can  be  solved  via  succes¬ 
sive  overrelaxation  (SOP)  if  the  right  hand  side  of 
the  equation  is  known.  For  this  study  the  simultane¬ 
ous  relaxation  method  with  overrelaxation  coefficient 
equal  to  1.0  was  used.  While  the  sequential  over¬ 
relaxation  technique  may  converge  quicker,  there  is  a 
limit  to  the  overreloxation  value  which  can  be  used  to 
maintain  computational  stability.  In  both  techniques, 
if  the  tolerane  for  convergence  is  set  to  low,  the 
number  of  iterations  becomes  large  and  computational 
instability  becomes  of  concern.  In  this  study  tests 
were  done  to  determine  where  computational  instability 
became  a  problem.  It  was  discovered  that  after  about 
100  iterations  utilizing  the  simultaneous  scheme 
results  began  to  degrade  due  to  computational 
problems.  Tolerence  levels  were  adjusted  accordingly 
to  insure  that  "acceptable  convergence"  was  reached 
prior  to  100  iterations.  The  convergence  criteria 
used  in  this  studv  were  found  to  be  lower  than  those 
of  other  studies  using  similar  techniques,  e.g. 


Wirsing  (1985),  Lin  e_t_  a_l .  (  198  5). 

Unfortunately,  in  the  right  hand  side  term  only  D 
is  known.  To  solve  (3.5.5),  the  anelastic  continuity 
equation  is  applied  to  the  vertically  integrated  hori¬ 
zontal  mass  divergence  so  that 

0  A  /  ^  _  A  \ 


p(&-&)d£= 


(3.5.7) 


Note  that  the  three  other  terms  on  the  right  hand 

_ A 

side  of  O.S.TJ.JWsFfcjf^P^SFc  are  all  zero  by  con¬ 
straint  and  assumption.  The  distribution  of  the  di¬ 
vergence  error  was  assumed  to  be  independent  of 
height.  With  the  aid  of  (3.5.6)  and  (3.5.7),  the 
field  of  Lagrange  multiplier,  A.  ,  could  be  determined 
assuming  the  natural  boundary  condition  (Sasaki, 
1970).  Then  horizontal  components  of  the  ’’true"  wind 
were  computed  from  (3.5.3)  and  (3.5.4).  Finally,  the 


vertical  velocity  at  any  height,  z^  ,  is  determined 
from  the  formula 

,  p.  ljUUi 


W  r  W  -  1  ifi - 

L  1  (3-5'8) 

The  derived  wind  field  not  only  satisfies  the  kine¬ 
matic  boundary  conditions  but  also  satisfies  the 
anelastic  continuity  equation.  This  internally  con¬ 
sistent  wind  field  is  then  used  to  recover  the  thermo¬ 


dynamic  variables. 


3.6  Computation  of  Thermodynamic  Variables 

Chong  e_t  al .  (  1980),  Kane  e_t  al .  (1931),  and  Gal- 
Chen  (1982)  showed  that  when  computing  perturbation 
quantities  such  as  perturbation  pressure,  density  and 
virtual  temperature,  it  is  absolutely  necessary  to  use 
a  storm- relative  coordinate  system.  When  dealing  with 
one  time  period,  this  essentially  becomes  a  correction 
for  storm  motion  and  scan  time  differences.  However, 
when  dealing  with  more  than  one  time  period,  salient 
features  such  as  downdraft  and  updraft  cores,  and 
reflectivity  patterns  can  be  used  to  orient  the  grid. 

In  mesoscala  meteorolgy,  it  is  important  to 
retain  the  internal  gravity  waves  propagating  through 
the  system,  while  at  the  same  time  filtering  out  the 
acoustic  or  sound  waves.  The  use  of  the  anelastic 
approximation  of  Ogura  and  Phillips  (1962)  in  pre¬ 
ference  to  Reynolds  form  of  the  momentum  equations 
achieves  this  result  and  at  the  same  time  simplifies 
the  computations.  The  anelastic  approximation  re¬ 
quires  the  hydrostatic-adiabatic  assumption  with 
density  varying  only  in  the  vertical. 

3.6.1  Computation  of  perturbation  pressure 

The  momentum  equations  in  a  moving  coordinate 
system,  see  Gal-Chen  (1978),  can  be  written  as 


(3.6.1) 


dn 


(3.6.2) 


(3.6.3) 


where 


f  =  -  +  (  '  W  0Ji<6  )  (3.6.4) 


.=  1,2,3 


Qj-~  ^V*Vv  V  ^  _  j^2  _TL  \A  S iV> p 

H- fV.Vw -( 3Jkv\ -  p^'CEji* f^y 


(3.6.5) 


aw  A,," 

t-  V** 


(3.6.6) 


( u  f  v , w )  are  velocity  components  in  the  x ,  y,  and  z 
direction,  respectively, 

^  is  the  basic  density  in  the  hydrostatic- 
adiabatic  atmosphere  which  is  approximately  equal  to 
the  environmental  density  (Brandes,  1984;  Lin  e_t  al .  , 
1985), 

/ 

■p  is  perturbation  density, 

/ 

P  is  perturbation  pressure, 

Tii-'tii.  ,  and^^  are  turbulent  momentum  fluxes, 
and  other  symbols  have  their  usual  meanings. 


The  turbulent  momentum  flux  terms,  , 


are  com¬ 


puted  using  the  formulae  by  Klemp  and  Wilhelmson 
(1978)  : 


and 


(3.6.7) 

(3.6.8) 

(3.6.9) 


from  Smagorinsky  (1963),  the  eddy 
ficient  is 


X 


PA  - 


t1  C \  2)uL 
\Jz  I 


viscosity  coef- 


(3.6.10) 


Next  we  parameterize  the  precipitation  drag  term,  , 
in  (3.6.6)  using  the  formulation  proposed  by  Douglas 


(1964)  which  relates  the  liquid  water  content  to  the 

Va. 


reflectivity  in  mm^/m^. 


where  a»1.82  and  b»2.4xl0* 


(3.6.11) 


Substituting  (3.6.7)  - 


(3.6.10)  into  (3.6.4)  -  (3.6.6)  we  finally  arrive  at 


<2  Pd  _  _  p\V-Vu_  +• 

3  %  a 


-s>,j 

J  5  ctr,  1  atj  +ri 


^  U,  SkJ la?,  ) 
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(3.6.12) 


where  Km  *»  eddv  momentum  mixing  coefficient,  C  =  0.21, 

Vs 

Cm  *  0.20,  and  Vi.  =  ( AX “as  suggested  by  Deardorff 
(  1975)  and  f  1  vsin#-wcos<#)  ,  f  2=*-f^2n.wcos#'  and 

f 3  =J^2Jlu cos#  are  the  Coriolos  forces. 

Once  the  turbulent  momentum  fluxes  have  been  com¬ 
puted,  the  horizontal  perturbation  pressure  can  be 
computed  by  using  Eqs.  (3.6.1)  and  (3.6.2).  Differen¬ 
tiating  (3.6.1)  with  respect  to  x  and  (3.6.2)  with 


respect  to  y  and  adding  yields 

if  +  sV _  if* a c. 

C>x^  "  C>x  'jy 


(3.6.15) 


This  Poisson  type  partial  differential  equation 
can  be  solved  using  sequential  relaxation  with  the 

retangular 


Neumann 


boundary  condition. 


For 
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boundaries,  we  have 


—  =  F 

;  ci-t  y  —  to. 

3y 
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Gal-Chen  (1978)  noted  that  a  unique  solution  to 
(3.6.15)  exists  only  if  the  horizontal  areal  mean  of 
P  ,  <P  >,  is  removed  from  P  of  (3.6.15).  The  devia¬ 
tion  perturbation  pressure  (P  *)  field  together  with 
the  wind  field  are  then  used  to  retrieve  deviation 
perturbation  density  and  temperature  fields. 

3.6.2  Perturbation  density 

Once  the  fields  of  deviation  perturbation 

pressure,  P^  ,  has  been  computed,  the  field  of  de¬ 
viation  perturbation  density,^'  (=  f7'  -  <p’>),  can  be 
obtained  from  Eq .  (3.6.3)  using  the  expression 

“  H"-  («-<*>)  (3.6.16, 

where  <?’>,  <  1 > ,  and  <H>  are  function  of  height 

only . 

3.6.3  Perturbation  virtual  temperature  t 

Once  the  fields  of  deviation  perturbation  j 

pressure,  ,  and  deviation  perturbation  density,  ^ ’ 


have  been  computed,  the  deviation  perturbation  virtual 
temperature  from  the  horizontal  mean, 


can  be  computed  using  the  perturbation  equation  of 


state 


-C  =  (fMa) 


(3.6.17) 


3.6.4 


fields 


Verification  of  the  derived  thermodynamic 


Since  the  method  of  determining  thermodynamic 
variables  presented  here  is  indirect,  verification  of 


the  derived  fields  is  essential 


Gal-Chen  and  Hane 


(1981)  suggested  a  verification  method  based  on 
momentum  checking  and  time  continuity.  The  momentum 
checking  method  involves  the  computation  of  a  dimen¬ 
sionless  quantity,  Er,  defined  as: 

jj(FV^)dxdy  (3'6-13> 

Er  will  be  computed  for  each  analysis  level.  Gal-Chen 
and  Hane  (1981)  noted  that  while  noise  produces  values 
of  Er  >  .5.  Only  values  of  Er  <  .5  can  be  interpreted 
scientifically.  However,  Hane  and  Ray  (1985)  pointed 
out  that  fields  which  produce  larger  values  of  Er, 
e .  g  .j  Er-s'l.O,  can  still  contain  usable  information. 
Er  as  a  measure  of  "goodness"  of  fit  of  the  derived 
perturbation  presssure  gradients  and  the  associated 
wind  field  should  be  looked  at  in  light  of  other 


checks  such  as  physical  linkage,  agreement  between 


derived  fields  and  vertical  velocity  or 
divergence/ convergence. 

The  time  continuity  check  involves  inspecting  the 
thermodynamic  fields  from  one  time  period  to  another. 
If  the  retrieved  thermodynamic  fields  maintain  good 
continuity  from  one  analysis  time  to  the  next,  the 
field  will  be  interpreted  as  reasonably  accurate. 


CHAPTER  4 


Error  Analysis 

Central  to  the  interpretation  of  the  results  of 
any  scientific  study  is  an  understanding  of  errors  in 
those  results  and  the  mechanisms  which  generated  them. 
Not  only  does  this  knowledge  aid  in  interpreting  the 
results  more  correctly  but  it  also  aids  in  future 
collection  and  processing  of  similar  data. 

Errors,  in  general,  are  random  or  non-random  in 
nature.  Those  which  are  random  are  related  to  obser¬ 
vational  techniques,  e.g.,  ground  clutter,  sxdelobe 
contamination,  sampling  problems,  temporal  errors  due 
to  sampling  differences  of  the  radars  involved,  sta¬ 
tistical  probability  considerations,  and  hardware 
processing  errors  among  other  things.  These  types  of 
errors  are  very  difficult  to  detect  and  almost  impos¬ 
sible  to  correct.  Non-random  errors  can  be  introduced 
from  a  number  of  sources,  e.g.,  errors  in  wind  field 
and  numerical  processing  techniques.  Numerical  pro¬ 
cessing  techniques  include  such  elements  as  choice  of 
boundary  conditions,  finite  difference  methods, 
direction  of  vertical  integration,  and  errors  in 
deduced  fields  due  to  errors  in  the  wind  field. 


4.1  Random  Errors 


Bohne  and  Srivastava  (1976)  and  Ray  et  al. 
(1978a)  examined  random  errors  associated  with  Doppler 
radar  wind  measurements.  Bohne  and  Srivastava  (1976), 
using  data  from  three  Doppler  radars,  deduced  that  the 
theoretical  accuracy  of  the  horizontal  wind  was 
approximately +1  m  s*  .  They  noted  that  the  accuracy  of 
the  terminal  velocity,  ,  is  about  +  2-3  m  s’’  for 
rain  and  +  1/2  m  s1  for  snow.  Ray  _et_  a_l.  (1978a) 
looked  at  the  differences  in  the  horizontal  wind 
fields  using  2,  3,  and  4  radars.  Distributions  of 
error  in  variances  for  two  radars  are  reproduced  in 
Fig.  4.1.  Note  that  for  the  two  radar  case,  a 
target  located  within  40  km  has  a  standard  deviation 
of  horizontal  wind  uncertainty  of  <  .2  m  s'  .  They 
note  that  increasing  the  number  of  radars  does  not 
appreciably  reduce  error.  Ray  e_t  a_l.  (1978a)  deduced 
that  errors  in  the  horizontal  wind  speeds  are  <3  m  s*. 
For  a  specific  storm  on  8  June  1974,  Doviak  et  al. 
(1976)  estimated  errors  in  wind  field  to  be  10  degrees 
in  direction  and  1  m  s^  in  speed.  Lin  and  Pasken 
(1982)  utilizing  a  formula  from  Lhermitte  and  Miller 
(1970)  arrived  at  similar  values.  Elmore  (1982)  con¬ 
ducted  an  error  analysis  suggested  by  Nelson  (1980) 
and  Nelson  and  Brown  (1982)  and  obtained  comparable 
results.  A  detailed  error  analysis  was  accomplished 
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on  results  of  this  study  and  is  discussed  in 
Chapter  5. 

Random  errors  in  velocity  measurements  are  asso¬ 
ciated  with  statistical  probability  and  are  difficult 
to  discuss  except  theoretically.  They  are  linked  to 
the  assumptions  made  about  the  drop  size  distribution 
within  the  scanned  volume,  the  composition  of  the 
hydrometers  within  the  scanned  volume,  and  sampling 
techniques.  One  example  involves  relating  the  terminal 
fall  speed  of  the  hydrometer  en  massed  within  the 
scanned  volume  to  the  refectivity.  Table  3.4  shows  a 
number  of  empirical  Z  -  V*  relationships  from  dif¬ 
ferent  authors  and  the  different  values  they  produce 
for  the  spectrum  of  reflectivity  values  expected  in 
this  study. 

Errors  associated  with  the  choice  of  Z  func¬ 
tion  are  generally  more  important  at  higher  elevation 
angles,  where  the  radial  component  of  the  wind  reveals 
a  larger  magnitude  of  the  vertical  wind.  In  this 
study,  elevation  angles  were  small  enough  that  this 
error  can  be  neglected. 

Some  of  the  errors  associated  with  hardware 
involve  dropped  bits,  incorrect  processing  of  the 
signal  due  to  low  signal  to  noise  ratio,  voltage 
fluctuations  within  the  radar  and  problems  similar  to 
that  involving  the  variance  measurements. 
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4.2  Non-random  errors 

Non-random  errors  can,  in  general,  be  divided 
into  two  subcategories: 

1)  Errors  associated  with  limiting  assumptions,  and 

2)  Errors  associated  with  numerical  techniques. 

4.2.1  Errors  associated  with  limiting  assumptions 

A  number  of  limiting  assumptions  were  made  in 
this  study.  Each  will  be  explained  in  the  following 
sub-sections . 

4. 2. 1.1  Quasi-steady  state  storm  assumption 

It  was  assumed  that  both  radars  were  scanning  the 
same  volume  at  the  same  time.  Great  effort  was  made 
during  JAWS  to  coordinate  the  radar  scanning  pro¬ 
cedures.  The  results  were  a  multiple  Doppler  radar 
data  set  which  has  temporal  differences  between  radars 
which  were  considered  negligible  for  this  study.  Of 
importance  is  the  problem  involving  reduction  of  the 
observation  time  needed  to  perform  a  complete  volume 
scan.  This  took  nearly  two  minutes.  The  solution  was 
to  correct  all  observations  to  a  "storm  time"  by 
correcting  the  location  of  the  observation  using  the 
time  correction  equations  in  Chapter  3.  Gal -Chen 
(1982)  suggested  that  this  correction  is  necessary  to 
reduce  errors  caused  by  storm  motion.  The  use  of 
(3.3.1)  and  (3.3.2)  involve  making 


two  assumptions 
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i 

|  about  the  state  of  the  storm.  The  first  assumption  is 

I 

j  that  for  the  entire  observation  period  the  motion  of 

i  the  storm,  as  a  whole,  was  constant.  Secondly,  not 

) 

!  only  was  it  constant,  but  it  was  moving  at  the  same 

» 

j  speed  at  all  levels.  Wilson  e_t_  al_.  (1984)  noted  for 

!  this  particular  storm  that  the  required  correction  was 

1 

negligible  for  a  grid  size  of  250  m.  Inspection  of 
j  the  results  from  Wilson  e_t  a_l .  (1984)  reveals  that 

!  changes  in  the  speed  of  this  storm  near  the  surface 

were  also  negligible.  The  maximum  time  differences 
j  between  the  "storm  time"  and  the  acquisition  time  for 

any  datum  point  is  less  than  80  seconds  which  suggests 
that  errors  in  location  of  a  datum  point  due  to  storm 
|  motion  is  less  than  one  grid  distance. 

|  The  quasi-steady  state  assumption  also  encom¬ 

passed  temporal  changes  in  the  structure  of  the  storm, 
j  Gal-Chen  (1978)  suggested  that  when  two  sequential 

scans  were  available,  the  sensed  velocity  at  a  given 
point,  after  having  been  corrected  spatially,  could  be 
j  corrected  temporally.  This  correction,  see  Eq.  (4.1), 

amounts  to  a  linear  interpolation  of  the  change  in 
velocity  using  the  first  and  second  scan  times  and  the 
storm  relative  time 

u  (  x,  y  ,  z  ,  t )  ;aJ0(u,r(  x  ,  y  ,  z  ,  tK)  +  u^  (  x  ,  y  ,  z  ,  t<  )  (4.1) 


> 
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The  time  derivative  can  be  evaluated  as 

!£L_T=  tj-  MT))/^  (4.2). 

This  procedure  requires  a  belief  that  when  all 
data  points  are  adjusted,  the  result  will  represent 
the  storm  as  it  appeared  at  the  "mean  observation 
time".  It  is  based  on  the  steady  state  assumption  for 
approximately  2  minutes.  This  technique  has  been 
employed  by  Lin  and  Pasken  (1982),  though  in  a  some¬ 
what  different  form,  and  Virsing  (1985)  in  a  similar 
f  orm . 


Lin  e_c  a^ .  (1985)  examined  the  effect  the  steady 
state  assumption  had  on  computations  of  deviation  per¬ 
turbation  pressure  and  deviation  perturbation  virtual 
temperature  and  concluded  that  the  errors  created  by 
making  the  steady  state  assumption  were  of  secondary 
importance.  They  also  noted  that  the  accuracy  of  the 
time  tendency  term  depends  on  the  "maturity"  of  the 
storm.  This  suggests  that  in  order  to  realistically 
include  the  time  tendency  term  s,  two  time  periods 
must  be  analyzed,  that  is,  parameterization  of  the 
time  tendency  term  is  questionable  without  knowing  the 
"age"  of  the  storm. 

Hane  e_t_  a_l.  (1981)  examined  the  effect  of  the 
steady  state  assumption  on  temperature  retrieval  using 
model  data.  Their  conclusion  was  that  while  the 


3 ' 

actual  results  deteriorated,  the  solution  was  judged 
useful  if  general  temperature  patterns  were  desired. 

4. 2. 1.2  Adiabatic-hydrostatic  assumption  on  density 
Table  3.5  shows  the  difference  for  densities  com¬ 
puted  from  the  Denver  sounding,  the  equation  of  state 
and  the  first  law  of  thermodynamics  making  the  adia¬ 
batic  and  hydrostatic  assumptions  necessary  for  the 
use  of  the  anelastic  continuity  equation.  This  dif¬ 
ference  has  little  effect  on  the  computation  of  the 
wind  field  due  to  differences  in  the  Foote  &  du  Toit 
correction,  see  Section  4. 2. 2. 2. 

4.2.2  Errors  associated  with  the  wind  field 

Some  of  the  errors  in  the  wind  field  have  been 
discussed  under  the  random  error  Section.  Non-random 
wind  errors  include  those  associated  with  specific 
formulations  of  the  terminal  velocity  equation,  cor¬ 
rections  of  the  terminal  velocity  using  the  Foote  and 
du  Toit  correction  scheme,  and  errors  associated  with 
subjective  analysis  of  the  horizontal  wind  field. 

4. 2. 2.1  Errors  due  to  terminal  velocity  assumptions 

In  this  study  we  have  chosen  to  use  the  terminal 
velocity-reflectivity  relationship  of  Martner  (1976). 
Table  3.4  shows  values  for  a  variety  of  reflectivity 
strengths  as  well  as  a  number  of  other  possible  equa¬ 
tions.  Examination  of  (3.11)  suggests  that  the  errors 


will  be  at  least  one  order  of  magnitude  smaller  than 
the  values  of  u,  and  v.  This  is  due  to  the  fact  that 


the  term  is  multiplied  by  z  -  ,  with  i  =  1  or  2, 
which  ranges  from  0-5  km,  while  x  -  ranges  from  0- 
30  km.  This  error  reaches  its  largest  magnitude  for 
the  higher  elevation  angles.  Table  3.4  suggests  that 
the  Martner  equation  provides  the  least  amount  of 
sensitivity  to  errors  which  may  be  associated  with  Z. 


4. 2. 2. 2  Errors  related  to  Foote  &  du  Toit  correction 
Foote  and  du  Toit  (1969)  suggested  two  forms  of  a 
correction  to  be  applied  to  the  computed  terminal 
velocities,  namely,  _ 

Y  ji 

V4  -  Voi  io  (i+aMZj(u-fXr-T.i)  {4  3) 


where 


and 


.  0.4- 


Y=  o.4iLo?|o(&J-0.4(Lc,!a(| 


2.5 


Vi,  (#) '  V, 


(4.4) 


where 


^  «  standard  surface  density  (Pe  -  1013.25  mb, 

T0  ■  293.16  K),  and  p  -  actual  density. 

This  correction  accounts  for  the  differences  in 
air  density  and  the  effect  of  this  difference  on 
aerodynamic  drag  of  the  droplets.  Figure  4.2  shows 
the  computed  vertical  velocities  using  (4.3)  and 
(4.4),  where  V  is  from  (3.4.2),  and  density,  P  ,  is 


TERMINAL  VELOCITY  (m/s) 


figure  4- .2.  Graph  of  terminal  velocity  for 

both  the  FULL ( . . . )  and  the  ABBREVIATED  ( _ ) 

Foote  4  duToit  correction  formula.  Reflectivity 
is  60d5t,  density  is  from  adiabatic-hydrostatic 
formulation.  Terminal  velocity  formula  is  from 
Hartner  (197.6) .  Difference  at  200  mb  is  1.08  m/s 
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from  (3.4.10).  While  these  differences  are  signifi¬ 
cant  when  related  to  the  actual  values  of  ,  in  light 
of  Section  4. 2. 1.2,  they  are  of  little  significance  in 
the  computation  of  the  u  and  v  components,  see 
(3.4.11a)  and  (3.4.11b). 

4. 2. 2. 3  Errors  associated  with  subjective  analysis 

There  is  no  exact  method  for  measuring  errors  in 
results  due  to  subjective  editing  of  the  data.  It 
suffices  to  say  that  for  this  study  wind  velocities  in 
the  storm  environment  were  modified  only  in  cases 
where  the  data  were  absolutely  required  to  perform  the 
study  . 

The  only  data  which  was  considered  for  modifica¬ 
tion  were  the  radial  winds  after  they  had  been  placed 
on  the  grid  and  prior  to  deriving  the  three-dimen¬ 
sional  wind  fields.  Among  the  considerations  prior  to 
modifying  a  given  piece  of  data  were;  1)  three- 
dimensional  continuity  of  the  modified  datum  within 
the  radial  wind  field,  2)  effects  that  the  modified 
data  had  on  the  derived  wind  fields,  and  3)  the  con¬ 
tinuity  of  the  derived  fields  with  and  without  the 
modifications. 


4.3  Errors  in  Derived  Thermodynamic  Fields 


Errors  in  the  pressure  field  come  basically  from 
two  sources.  Errors  in  the  wind  field  will  produce 


errors  in  the  retrieved  pressure  field.  Pasken  (1981) 
used  model  output  from  the  Lin  and  Chang  (1977)  storm 


model  to  examine  this  type  of  error.  He  found  that 
10%  error  in  the  wind  field  produced  approximately  a 
10%  error  in  the  retrieved  pressure  field. 

The  actual  output  from  the  pressure  retrieval 
program  is  the  perturbation  pressure  field.  This 
solution  is  not  unique.  Removing  the  mean  perturba- 


tion  pressure  from  the 

solution 

field 

leaves 

the 

deviation  perturbation 

pressure 

field, 

which 

is 

unique.  If  the  mean  perturbation  pressure  which  is 
removed  is  in  error,  this  will  result  in  an  error  in 
the  deviation  perturbation  pressure  fields.  Knowledge 
of  the  mean  perturbation  pressure  is  not  available, 
however,  computations  of  the  mean  and  standard  devia¬ 
tions  of  the  resultant  pressure  fields  were  done  level 
by  level.  If  the  area  of  computation  of  the  pressure 
field  was  large  enough,  it  was  assumed  that  the  mean 
of  the  perturbation  pressure  would  be  approximately 
zero  (i.e.,  <p’>^0,  Gal-Chen,  1978).  Examination 

of  the  deviation  perturbation  means,  not  shown,  reveal 
that  these  values  are  zero.  It  suggests  that  the 
deviation  perturbation  pressure  pattern  is  realistic 
and  probably  free  of  error  of  the  second  type.  Lin  e_t 
a  1  .  (  1985)  investigated  the  sensitivity  of  pressure 

retrieval  to  exclusion  of  the  local  time  tendency. 


They  discovered  that  the  inclusion  of  the  local  time 
derivative  terms  in  the  dynamic  equations  did  not 
improve  the  accuracy  of  pressure  recovery. 

Errors  in  temperature  retrieval  depend  heavily  on 
the  input  wind  fields  as  well  as  the  retrieved  pres¬ 
sure  and  density  fields.  The  use  of  friction  in  the 
computations  of  both  pressure  and  temperature  ulti¬ 
mately  involve  the  second  and  third  derivative  of  the 
wind  field.  Gal-Chen  and  Kropfli  (1984)  investigated 
the  effect  of  friction  in  the  planetary  boundary  layer 
on  the  retrieval  of  pressure  and  potential  temperature 
fields.  They  found  that  the  inclusion  of  friction  did 
not  significantly  reduce  relative  errors  (Er)  in 
pressure  recovery  and  suggested  that  the  pressure 
calculation  seems  to  be  relatively  insensitive  to 
which  frictional  parameterization  scheme  is  used. 

Using  model-generated  data  as  input,  Hane  ejt  al . 
(1981)  inserted  random  errors  of  +  0.5  m  s'  into  the 
wind  field.  Their  findings  indicate  that  retrieved 
perturbation  virtual  temperature  is  sensitive  to 
errors  in  the  wind  field.  They  found  that  errors 
introduced  in  the  form  of  "white  noise",  implying 
wavelength  independence,  produced  an  average  error  in 
the  associated  temperature  field  of  about  1.0  C,  on 
the  average,  see  their  Fig.  12.  They  further  sug¬ 
gested  that  if  the  wavelengths  at  which  the  noise  was 


CHAPTER  5 


Sensitivity  Tests  of  Thermodynamic  Retrieval 

A  series  of  experiments  were  conducted  in  an 
effort  to  determine  the  sensitivity  of  the  thermo¬ 
dynamic  retrieval  method.  Many  investigators  have 
looked  at  whether  the  inclusion  of  a  correction  scheme 
which  brings  the  continuity  equation  into  balance 
after  correcting  the  vertical  velocities,  aids  in 
producing  "core  correct"  results,  e.g.,  Ray  e_t  a  1  . 
(I960),  Ziegler  e_t  a_l  .  (1983),  Gal-Chen  and  Kropfli 
(1984)  and  Lin  e_t  a_l.  (  1985).  In  addition,  Lin  e_t  a  1 . 
(1985)  and  Wirsing  (1955)  looked  at  the  resultant 
thermodynamic  fields  when  no  correction  was  imposed. 
In  both  these  cases,  downward  integration  was  used 
from  the  storm  top  and  the  storm  was  considerably 
stronger  than  the  storm  investigated  in  the  present 
study.  The  overwhelming  conclusions  from  these  and 
other  investigations  are  that  some  type  of  correction 
of  the  3D  wind  field  is  necessary  in  order  to  satisfy 
mass  continuity  within  the  domain  of  interest.  This 
correction,  based  on  the  variational  concept,  is  par¬ 


ticularly  needed  if  observational  errors  are  signifi- 


cant.  To  this  end,  sensitivity  of  the  recovery 
technique  to  inclusion  of  a  correction  scheme  was  not 
tested.  All  results  presented  herein  include  the 
calculus  of  variation  scheme  (Sasaki,  1970)  as  the 
correction  scheme  unless  otherwise  noted. 

In  this  study,  we  attempted  to  further  examine 
the  impact  of  random  errors  in  the  wind  field  on  the 
calculation  of  vertical  velocity,  pressure  and  temper¬ 
ature  fields  for  a  non-severe  convective  storm.  The 
following  cases  were  considered: 

1)  Inclusion  of  two  types  of  frictional  parameteriza- 
tions  in  the  momentum  equations. 

2)  Smoothing  of  the  radial  wind  field  by  increasing 
the  scan  radius  in  the  3 a r nes'  objective  analysis 
scheme . 

3)  Sensitivity  of  the  retrieval  process  to  randomly 
introduced  errors  in  the  initial  wind  field. 

4)  Sensitivity  to  changes  in  liquid  water  content. 
Table  5.1  describes  the  identification  of  the  cases 
and  the  elements  included  or  excluded  from  each  of  the 
experiments. 

Gal-Chen  (1978)  suggested  a  measure  of  "goodness" 
of  the  recovery  technique.  The  momentum  check  men¬ 
tioned  in  Chapter  3  provides  a  measure  of  "relative 
error"  between  the  derivatives  of  the  recovered  pres¬ 
sure  field  and  the  momentum  equations,  see  Eq  . 


(3.6.13).  Small  values  of  ' 
better  agreement  between  tl 
momentum.  However,  this  coi 
but  not  sufficient.  For  th: 
to  examine  the  physical  pi; 
tinuity  of  the  retrieved  fi< 
Hane,  1981).  According  to  H; 
values  still  allow  for  extr; 
tion  from  the  recovered  f: 
employed  in  this  study 
included  or  excluded  term  w< 


wavelength  of  a  minimum  resolvable  wave  as  discussed 
in  Appendix  A.  Therefore  it  plays  a  role  somewhat 
similar  to  smoothing  discussed  in  studies  by  Gal-Chen 
and  Kropfli  (1934),  Lin  e_t  a_l  .  (1985)  Wirsing  (1935), 
etc.  The  studies  by  Wirsing  (  1985)  and  Lin  ej^  a  1  . 
(1985)  discovered  that  smoothing  substantially  reduced 
the  momentum  check  values.  Unlike  an  arbitrary 
smoothing  scheme  (Shuman,  1959)  used  in  the  above  two 
studies,  the  scheme  employed  in  this  study  allows  us 
to  explicitly  know  what  is  the  minimum  resolvable 
wavelength  in  the  derived  fields.  More  detailed 
discussions  follow. 

Modification  of  the  liquid  water  content  produced 
deceiving  results.  When  viewed  from  the  perspective 
of  area  averages,  the  recovery  process  seemed  to  be 
relatively  insensitive.  However,  sensitivity  appears 
to  be  directly  related  to  reflectivity  and  can  become 
significant  where  reflectivities  are  strong. 

5.1  Frictional  Parameterization 

Two  frictional  parameter izations  were  used  in 
this  study.  The  first  was  based  on  the  Austausch 
hypothesis  with  constant  eddy  viscosity  coefficients, 
labeled  A  in  the  text.  Both  horizontal  and  vertical 
coefficients  were  set  to  be  100  m^s  '  .  The  second 
approach  was  suggested  by  K  letup  and  Wilhelm son 
(1973),  labeled  KW  in  the  text.  The  governing  equa- 


tions  for  this  parameterization  can  be  found  in 
Section  3.6.  Table  5.2  shows  the  momentum  check  (2r) 
results  for  the  various  cases.  Examination  of  VUN 
versus  VUKW  reveals  an  average  reduction  of  the  momen¬ 
tum  check  values  of  less  than  5%.  However,  much 
improvement  is  found  when  case  V S H  is  compared  to  case 
VSXV.  It  appears  that  for  unsmoothed  data  frictional 
parameterization  has  less  effect  than  for  smoothed 
data  in  pressure  retrieval.  Table  5.2  further  shows 
that  when  comparing  either  frictional  case  (VS A  or 
VSXW)  with  the  smoothed  case  without  friction  (VS!.'), 
inclusion  of  friction  is  significant  in  the  lower 
layers  but  not  above  4.5  km.  The  larger  momentum 
check  value  near  the  surface  is  the  result  of 
numerical  problems  with  a  one  sided  derivative  com¬ 
pounded  by  lower  boundary  condition  that  vertical 
velocity  vanishes  at  the  surface.  Inspection  of  the 
reflectivity  pattern  and  the  sounding  reveal  that  the 
cloud  base  of  the  parent  cloud  was  near  4.0  km  AGL. 

As  previously  stated  the  vertical  and  horizontal 
eddy  viscosity  coefficients  were  both  set  to  be 
100  m2  s-1  The  similarity  between  VSK !■/  and  V3A  is 
evident.  This  would  indicate  that  while  inclusion  of 
some  friction  parameterization  is  needed,  it  does  not 
seem  to  matter  which  scheme  is  used.  This  result 
agrees  with  Hane  et  al.  (1931).  They  concluded  that 
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exclusion  of  the  frictional  parameterization  would 
degrade  results  but  that  the  solution  should  be  of 


value.  These  results,  however,  do  not  agree  with 
other  investigations,  e.g.,  Gal-Chen  and  Xropfli 
(1984)  or  Lin  e_t  al .  (  1985).  The  study  by  Lin  ej:  a  1 . 

(1935)  found  that  friction  did  not  contribute  signi¬ 
ficantly  to  thermodynamic  retrievals  at  any  level 
inside  the  severe  storm  embedded  within  a  squall  line 
system..  The  lowest  level  analyzed  in  that  study  was 
1.5  km  AGL  and  it  appeared  to  be  above  the  base  of  the 
parent  cloud.  This  would  tend  to  suggest  that 
frictional  contributions  nay  be  most  important  in  the 
area  below  the  cloud  and  within  first  1  or  2  km  of  the 
surface.  On  the  other  hand,  Gal-Chen  and  Xropfli 
(1984)  studied  the  dynamics  of  the  planetary  boundary 
layer  for  a  clear  sky.  They  found  no  significant 
reduction  in  the  relative  errors  (Hr)  for  the 
retrieval  of  pressure  and  temperature  perturbations 
when  the  subgrid  scale  turbulence  parameterization  was 
included  in  the  dynamic  equations. 

Unlike  the  above  two  studies,  the  present  study 
is  focused  on  the  structure  of  a  nonsevere  convective 
storm  in  Colorado.  This  type  of  storm  usually  has  a 
deep  subcloud  layer  up  to  4  km  AGL.  During  the  occur¬ 
rence  of  a  microburst  in  the  lower  troposphere,  strong 
3D  wind  shears  develop  causing  severe  turbulence.  As 


a  result,  a  contribution  due  to  convergence/divergence 
of  momentum  transfer  by  eddies  to  the  momentum  change 
becomes  significant.  This  is  evident  in  the  Er  values 
presented  in  Table  5.2.  Note  that  at  most  levels  the 
Austausch  formulation  produces  slightly  better  results 
than  the  K 1 e m p  and  Wilhelmson  (1978)  scheme.  The 
standard  deviations  of  pressure  ( S  D  p  * )  are  higher, 
especially  above  2  km,  when  the  Klenp  and  Wilhelmson 
(197S)  scheme  is  used,  while  the  standard  deviations 
of  temperature  ( S DT ' )  are  higher  when  the  Austausch 
formulation  is  used,  see  Tables  5.3a  and  5.3b.  It 
appears  that  while  the  contributions  from  friction  to 
retrieved  pressure  nay  be  inversely  proportional  to 
TT,  contributions  to  temperature  appear  to  be  propor¬ 
tional  to  X  , 

In  order  to  learn  more  about  the  effect  of 
friction  on  the  dynamics  of  a  micro  burst-producing 
nonsevere  storm,  spatial  variation  of  eddy  viscosity 
coefficients  (Km)  and  Reynolds  stresses  (Z )  using  the 
frictional  parameterization  of  Klenp  and  Wilhelmson 
(1978)  were  investigated.  Results  are  presented  in 
Figs.  5.1  through  5.5  and  Table  5.4. 

Table  5.4a  shows  the  means  and  standard 
deviations  of  Reynolds  stresses  (X  )  created  using  the 
Klenp  and  Wilhelmson  (1973)  scheme.  Table  5.4b  shows 
similar  components  of  Reynolds  stresses  (X")  from  the 


Table  5.3a 


Values  of  the  standard  deviation  (SD)  for  the  two 
retrieved  thermodynamic  variables  at  a  given  level  for 
1647  MDT.  Cases  considered  include  unsmoothed  no 
friction  (VUN),  Klemp  and  Wilhelmson  frictional 
parameterization  (VUKW),  and  Austausch  frictional 
parameterization  (VUA)  with  Km  =  100  The  domain 
being  considered  is  10  km  by  10  km  centered  on  the 
microburst  and  containing  the  entire  storm.  Units  are 
in  Pascals  (Pa)  and  for  pressure  and  temperature, 
respectively  . 


(VUX) 


(VUKW) 


(VUA) 
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Table  5.4a 


Values  of  the  mean  and  standard  deviation  (SD)  for 
horizontal  Reynolds  stresses;  't' xx,  'fxy  and  'ZSyy  and 
vertical  Reynolds  stresses;  *£Txz,  /T2yz,  and  zz  , 
computed  from  case  VSKW  with  the  Klemp  and  Wilhelmson 
frictional  formulation  at  a  given  level  for  1647  MDT . 


Horizontal 


i  xx 


L  xy 


r 


yy 


Z(ka) 

mean 

SD 

mean 

SD 

mean 

S  D_ 

N/m* 

N  /  m"2. 

H/n1 

N/o2, 

N/ma 

N/m** 

0.25 

-.013 

.  106 

-.001 

.042 

-.024 

.152 

0.50 

-.015 

.094 

-.006 

.041 

-.017 

.099 

0.75 

-.007 

.048 

- .  006 

.027 

-  .019 

.070 

1.0 

-.007 

.  03S 

-.006 

.032 

-.040 

.132 

1.5 

-.024 

.078 

-.004 

.053 

-.033 

.115 

2.0 

-.033 

.107 

-.004 

.059 

-.029 

.104 

2.5 

-.042 

.112 

-.011 

.  056 

-.039 

.139 

3.0 

-.054 

.124 

-.013 

.042 

-.061 

.187 

3.5 

-.067 

.139 

-.010 

.061 

-.045 

.147 

4.0 

-.083 

.150 

-.009 

.054 

-.040 

.122 

4.5 

-.082 

.157 

-.008 

.047 

-.035 

.135 

5.0 

-.056 

.134 

-.001 

.042 

-.026 

.081 

5.5 

-.039 

.124 

.006 

.048 

-.037 

.173 

6.0 

-.032 

.118 

.008 

.067 

-.033 

.175 

6.5 

-.026 

.094 

.006 

.060 

-.022 

.107 

7.0 

-.020 

.099 

.005 

.  065 

-.008 

.066 

7.5 

-.008 

.077 

.003 

.053 

-.003 

.051 

8.0 

-.006 

.071 

-.002 

.042 

-.004 

.050 

8.5 

-.006 

.067 

-.002 

.029 

-.010 

.062 

QO 


Table  5.4a 
continued 


As  in  Table  5.4a  but  for  case  VSA  with  constant 
viscosity  coefficient  (Km  =  100  m2s-'  ). 

Horizontal 

Txx  Xu 

Z(km)  mean  SD_  mean  SD  mean 


Table  5.4b 
continued 


Vertical 


l  xz 


mean  SD  mean  SD 


mean 


Austausch  scheme.  Figures  5.1  and  5.2  show  the  mean 
of  layer  averages  of  horizontal  Reynolds  stresses  for 
case  VS \'J  and  case  VS  A,  respectively.  A  comparison  of 
these  components  reveals  that  while  they  are  similar 
in  appearance,  for  the  Austausch  formulation,  the 
means  of  the  T  values  are  more  positive  and  tend  to 
be  slightly  smaller  than  those  generated  by  the  Klemp 
and  V.’ ilhelmson  (1973)  formulation.  Figure  5.3  depicts 
the  reason  for  this  phenomenon. 

Figure  5.3  shows  the  computed  x  x  for  various 
values  of  L'C  i-L  X .  The  horizontal  axis  represents  the 
value  L T  and  assumes  =  500  m.  For  a  span  of  £.0 
from  -10  to  +10  representing AU/AX  of  -.02  to  +.02., 
the  values  produced  by  the  Austausch  formulation  are 
seen  to  be  smaller  than  those  from  the  Klemp  and 
W  ilhelmson  formulation.  When  the  XI.  emp  and  V  ilhelmson 
parameterization  is  used  to  compute  the  Reynolds 
stress,  the  values  are  asymmetrical  about  0.  That  is, 
a  negative  value  of  A’j  /A  X  produces  larger  (in  the 
absolute  sense)  values  of  than  the  corresponding 
positive  value  of  AU/A-X.  This  accounts  for  the  nega¬ 
tive  shift  of  the  Reynold  Stress  values  when  the  Klemp 
and  V ilhelmson  parameterization  is  used.  Similar 
results,  not  shown,  occur  for  the  vertical  Reynolds 
stress  components.  h'ote  that  the  computations  of  the 


KM  <AGL> 


ecsr 


5.1  Horizontal 
or  case  VSIC’.v. 


( /  n2  ) 


Reynolds  stresses  TSx x  ,  'L-'t y  and 
Units  are  Newtons  per  square 


Fij*  5.3  Range  of  u xx  for  given  U— component  dif 
f erences  for  klenp  and  '/  i  1  h  e  1  m  s  o  n  frictional  para 
me t er i za t io n  and  Austausch  frictional  para 
meter izat ion . 


absolute  value  of  L,  .  Since  these  means  are  so  close 
to  zero,  the  standard  deviations  are  very  close  to  the 
means  of  the  absolute  values  of  ^  •  Actual  values  of 
in  the  vicinity  of  the  microburst  center,  not 
shown,  ranged  in  value  from  .3  to  -.3  for  the  VSKW 
case  • 

Figures  5.4  and  5.5  show  the  mean  horizontal 
(Kmxx,  Kmxy  and  Kmyy)  and  vertical  (Kmxz,  Kmyz  and 
Kmzz)  eddy  viscosity  coefficients,  respectively,  for 
case  VSKW.  The  larger  values  of  Kmyz  found  in  the 
raid-troposphere  are  related  to  the  v  component  jet 
which  was  previously  identified.  These  findings 
suggest  that: 

1)  When  compared  to  the  frictional  formulation  based 
on  the  Austausch  hypothesis,  the  Klemp  and  Wilhelmson 
formulation  produces  lower  estimates  of  the  frictional 
contribution  in  the  pressure  retrieval  process  by 
underestimating  the  eddy  viscosity  coefficient. 

2)  When  compared  to  the  Austausch  formula,  Klemp  and 
Wilhelmson  frictional  scheme  tends  to  smooth  the 
recovered  temperature  and  intensify  pressure 
deviations . 

5.2  Smoothing  of  Radial  Winds 

Smoothing  in  this  study  was  accomplished  by 
increasing  the  scan  radius  in  the  Barnes  gridding 
procedure  from  1.25  to  1.75  km.  Prior  to  the  gridding 
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<KJm>  CM  VS  3 

VOTICAL  EDDY  VISCOSITY  COEF.  <VSKV> _ 

Fig.  5.4  Vertical  eddy  viscosity  coefficients 
computed  from  basic  case  (VStiW)  horizontal  winds 
using  formula  of  Sraagorinslcy  (1963).  Units  are 
^s''. 
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process,  the  data  were  sliced  into  slabs  which,  except 


for  the  levels  below  1.0  km,  were  mutually  exclusive, 
that  is  data  from  one  layer  were  not  used  in  the  grid- 
ding  process  of  another  layer.  Thus  increasing  the 
scan  radius  essentially  resulted  in  a  horizontal 
smoothing  rather  than  a  full  three  dimensional 
smoothing  process.  Table  5.2  shows  the  results  of  the 
momentum  check  values  for  the  unsmoothed  no  friction 
case  (VUN)  and  the  smoothed  no  friction  case  (VSN). 
Note  that  the  percentage  decreases  in  Er  are 
particularly  significant  in  the  layer  from  the  surface 
to  4  km.  This  is  approximately  the  height  of  the 
parent  cloud  base.  Table  5.5  shows  the  means  and 
standard  deviations  of  the  three-dimensional  wind 
components.  As  expected  the  means  and  standard 
deviations  for  the  smoothed  case  ( VSN  )  are  less  than 
for  the  unsmoot'ned  case  (VUN').  Notice  that  the 
apparent  wind  maximum  in  the  v  component  occurs  in  the 
layers  between  4  to  6  km.  Due  to  the  mutual  exclu¬ 
sivity  of  data  used  in  the  level  gridding  procedure, 
existence  of  this  feature  in  the  vertical  indicates 
that  this  is  an  actual  feature  rather  than 
contamination  of  a  number  of  grid  points  in  the  hori¬ 
zontal  by  an  error  at  a  single  data  point. 

As  noted  earlier,  Table  5.3  shows  the  standard 
deviations  for  the  pressure  and  temperature  variables 
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retrieved.  Results  show  th  t  smoothing  (VS.'I)  reduces 
the  standard  deviations  when  compared  to  the 
unsmoothed  case  (  V  U  N  )  .  This  is  to  be  expected  since 
to  a  great  extent  these  retrieved  variables  depend  on 
the  wind.  Further,  in  both  smoothed  and  unsmoothed 
results  there  appears  to  be  a  maximum  located  within 
the  jet  level.  Examination  of  the  individual  terms  in 
the  momentum  equations  reveals  that  most  of  the 
variability  comes  from  the  S  v/^z  term. 

5.3  Inclusion  of  Random  Error 

Tha  actual  amount  of  errors  in  the  measured  wind 
field  comes  from  several  sources,  see  Chapter  4  for 
details.  however,  it  has  been  estimated  that  these 
errors  are  on  the  order  of  IOC  of  the  measured  wind, 
e.g.,  Waldteufl  (1976),  Doviak  and  Zrnic  (19S4),  etc. 
Since  no  independent  measurements  were  made  of  this 
storm,  it  is  impossible  to  determine  actual  errors 
present  in  the  derived  winds  or  retrieved  variables. 
It  is,  however,  possible  to  determine  sensitivity  of 
the  random  error  inherent  in  the  wind  field  to  the 
retrieval  of  pressure  and  temperature  perturbations. 
We  have  chosen  the  case  V  S  X W  to  conduct  the 
sensitivity  experiment.  Two  tests  were  performed  . 
In  each  test  random  error  amounting  to  a  maximum  of 


10°)  (VSKWE1)  or  20"  (VS5CWE2)  of  the  gridded  radial 
wind  for  each  radar  was  inserted  into  the  wind  field. 


The  data  were  then  processed  in  a  manner  similar  to 
case  VST'.'  discussed  previously. 

Figures  5.6  and  5.7  show  the  retrieved  perturba¬ 
tion  pressure  and  temperature  fields,  respectively,  at 
0.25  km  for  case  V  S  I\  W  E  1  .  The  pressure  and  temperature 
patterns  for  case  V S T E 1  (Figs.  5.6  and  5.7)  are  quite 
similar  to  those  for  the  base  case  V S \ W  (Figs.  6.14 
and  6.19).  This  indicates  that  the  contaminated  wind 
field  did  not  greatly  distort  the  retrieved  pressure 
and  temperature  fields.  Some  useful  information  still 
can  be  extracted  from  t.e  retrieved  fields.  This 
finding  is  encouraging  since  it  demonstrates  that 
thermodynamic  retrieval  is  not  quite  sensitive  to  a 
typical  error  of  up  to  10"  in  the  derived  wind  field. 
For  this  reason,  we  believe  that  the  retrieved  fields 
of  deviation  perturbation  pressure  and  temperature  are 
quite  reliable.  Table  5.6  shows  the  root  mean 
square  errors  (RUSEs)  and  the  standard  deviations  of 
the  differences  between  the  derived  pressure  and 
derived  temperature  fields  which  are  free  of  induced 
error  ( V  S  K  W )  and  those  with  induced  error  (VST'./  El)  and 
(VSTVE2).  Notice  that  magnitudes  of  RMSE  for  the  test 
with  10.1  error  (  case  VSTWEl)  are  relatively  small  in 
comparison  with  those  of  2  0  Z  error  (case  V  S  X E  2  )  . 
This  rasult  indicates  that  the  thermodynamic  retrieval 
technique  is  viable  if  the  derived  wind  field  is 
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1  .■  . . 


T»  . 


1  17 

within  10  to  157  accuracy  of  the  "true"  wind  field. 

The  sensitivity  of  the  retrieval  process  can  be 
measured  by  comparing  the  standard  deviations  of  the 
perturbation  pressure  and  temperature  from  case  V  S  E  W  , 
presented  in  Table  5.3b,  to  the  root  mean  square 


errors 

shown  in 

Table  5.6 

.  Results 

show 

that  errors 

created 

in  the 

retrieved 

pressure 

and 

temperature 

fields 

are  of 

the  same 

order  of 

magnitude  as  the 

induced 

errors, 

i  .  e  .  ,  for 

the  10. 7  induced 

error  case, 

there  appears  to  be  about  107  differences  in  the 
retrieved  field.  This  case  represents  the  magnitude 
of  the  error  expected  in  the  wind  field  synthesized 
from  Doppler  data.  For  case  VSK'.v'EE,  a  207  induced 
error  seems  to  produce  errors  of  about  207  to  407. 
N'ote  that  inherent  in  this  comparison  is  the  fact  that 
the  "base"  pressure  and  temperature  fields  are  not 
absolutely  free  of  error  themselves.  Errors  exist  in 
the  "base"  wind  field  used  to  produce  the  results  of 
case  VS KW.  Thus  any  suggestion  that  there  is  a  linear 
relationship  between  the  magnitude  of  errors  in  the 
derived  wind  field  and  the  retrieved  thermodynamic 
fields  should  be  thought  of  as  relative  only  and  not 
absolute.  These  results  are  similar  to  those  reported 
by  Pasken  (1981)  where  similar  comparisons  were  made 
using  model  data.  Pasken  (1931)  concluded  that  there 
seems  to  be  a  direct  relationship  between  the  amount 


of  inserted  error  and  the  resultant  error  on  the 
retrieved  fields. 


5.4  Sensitivity  to  Liquid  Water  Loading 

Sensitivity  to  liquid  water  content  changes  were 
investigated  by  Hane  e_t  al_.  (1981).  They  concluded 
that  doubling  the  liquid  water  content  did  not 
contribute  significantly  to  degradation  of  the 
retrieved  variables.  A  similar  test  was  performed  in 
this  study.  The  procedure  was  to  double  the  liquid 
water  content  (VSKWEL)  in  the  precipitation  drag 
computation  in  the  temperature  retrieval  process. 
Since  this  process  did  not  affect  the  pressure 
computation,  comparison  of  the  momentum  checks  cannot 
be  made  since  both  (VSKW)  and  (VSKWEL)  used  the  same 
derived  pressure  fields. 

Column  5  in  Table  5.6  shows  the  RMSE  values  of 
temperature  retrieval  between  cases  VSKUEL  and  VSKW. 
When  these  values  are  compared  with  the  standard 
deviations  of  temperature  for  case  VSKW,  the  RMSE 
values  are  much  smaller  in  magnitude.  It  is  apparent 
that  this  scheme  is  somewhat  insensitive  to  changes 
in  the  liquid  water  content.  Figure  5.3  shows  the 
temperature  differences,  T-Te,  for  0.5  km  level. 


Results  show  that  errors  in  liquid  water  content  tend 
to  be  related  to  the  reflectivity  pattern.  Although 
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the  overall  average  changes  indicate  small  changes, 
Fig.  5.3  shows  that  the  errors  can  be  very  significant 
in  areas  with  high  reflectivity.  In  addition,  the 
same  test  was  performed  with  the  precipitation  drag 
term  removed  from  the  computation  (not  shown).  The 
results  were  comparable  to  those  when  the  liquid  water 
content  was  doubled.  The  inclusion  of  the  precipita¬ 
tion  drag  term  in  the  computation  of  deviation  pertur¬ 
bation  (virtual)  temperature  seems  to  be  warranted  by 
these  results.  This  would  be  especially  true  if  the 
microburst  being  investigated  were  wet,  where  precipi¬ 
tation  drag  term  plays  an  important  role  in  the 
temperature  computation. 

5.5  "Best"  Retrieval  Method 

From  the  results  presented  in  Tables  5.2  through 
5.6,  Figs.  5.1  through  5.8,  and  using  the  momentum 
check  as  a  measure  of  goodness,  the  smoothed  case 
including  friction  (VSA  or  VSKV)  produces  the  better 
results.  Since  the  Klenp  and  Wilhelmson  (1973) 
formulation  appears  to  be  physically  more  acceptable, 
the  results  of  this  case  will  be  discussed  in 


Chapter  6. 


CHAPTER  6 


Discussion  of  Results 

In  this  study  ten  experiments  were  conducted  to 
ascertain  the  sensitivity  of  the  retrieval  process  to 
various  modifications  of  a  basic  process,  see  Table 
5.1.  Findings  show  the  following: 

1)  Correction  of  the  horizontal  wind  to  balance 
the  anelastic  continuity  equation  is  needed. 

2)  Smoothing  in  the  form  of  increased  scan 
radius  used  during  the  Barnes'  gridding  process 
produces  better  results. 

3)  Inclusion  of  a  frictional  parameterization 
decreases  momentum  check  values. 

4)  Inclusion  of  a  precipitation  drag  term 
improves  results. 

All  of  these  processes  were  incorporated  into  case 
VSKW  and  the  results  are  discussed  below. 

6.1  Three-dimensional  Flow  Pattern 

6.1.1  Horizontal  velocity 

Figures  6.1  and  6.2  show  the  pattern  of  u  and  v 
velocity  components , respectively ,  at  0.25  km  for  1647 
MDT .  For  comparison,  storm  relative  wind  vectors  are 

plotted.  All  heights  are  above  ground  level 
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also 
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(AGL).  Grid  distances  are  in  kn  with  the  origin 
located  on  the  mi c r o bu r s t .  The  center  of  the  micro- 
burst  not  shown  in  the  figures  is  located  at  12.5  km 
west  and  13.5  km  south  of  the  CP-2  radar  (see  Fig. 
3.4).  Line  AB  in  Fig.  6.1  denotes  the  location  of 
vertical  cross  section  presented  later.  This  vertical 
cross  section  is  in  line  with  the  environmental  flow 

vector.  Contours  for  Figs.  6.1  and  6.2  are  every 

- 1 

2  n  s  .  The  center  of  the  microburst  can  be  located 
by  superimposing  the  zero  contour  from  Fig.  6.1  onto 
Fig.  5.2  (dashed  line).  The  intersection  of  the  two 
zero  lines  represents  the  center  of  the  micro  burst  at 
0.25  km.  Note  that  the  absolute  wind  difference  in 
the  u  component  across  the  microburst  is  about 
10  m  s  and  slightly  larger  than  that  for  the  v  com¬ 
ponent.  The  criterion  for  a  downburst  to  be  a  micro¬ 
burst  is  met  by  both  components  of  the  horizontal 
wind,  see  Fujita  (1978)  and  McCarthy  e_t  a_l .  (1933). 

The  shallow  extent  of  the  horizontal  wind  shear 
can  be  seen  by  comparing  Fig.  6.1  at  0.25  km  to  Fig. 
6.3  representing  the  same  diagram  but  for  1.0  km  AGL. 
At  1.0  km  the  u  component  of  wind  is  negative 
(southerly)  across  the  entire  grid.  Absence  of  the 
outflow  feature  especially  in  the  u  component  of  hori¬ 
zontal  wind  suggests  that  the  horizontal  outflow  from 
this  microburst  is  confined  in  a  shallow  layer  from 
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the  surface  to  1  km  AGL.  This  is  consistent  with  many 
of  the  examples  cited  by  Fujita  (1535). 


Figure  6.4  shows  the  reflectivity  pattern  at  0.25 
km  superimposed  on  the  horizontal  wind  field.  The 
center  of  the  microburst  is  marked  by  "M".  Contours 
are  10  d3z  beginning  with  20  dBz.  The  high  reflec¬ 
tivity  within  0.25  kn  of  the  surface  suggests  that 
this  is  a  wet  microburst.  Wilson  e_t  a_l_ .  (  1 9 S 4 )  dis¬ 
played  the  0.05  km  pattern  of  reflectivity  and  wind 
for  1541,  1643,  1646  and  1643  MDT.  They  showed  a 
similar  high  reflectivity  pattern  50  n  from  the  sur¬ 
face  for  the  entire  period.  This  lends  credence  to 
the  suggestion  that  at  1647  MDT  this  storm  was  in  its 
mature  stage  rather  than  developing  and  was  a  wet 
rather  than  dry  microburst.  It  is  unfortunate  that 
this  storm  occurred  outside  the  PAM  sampling  region. 
Thus  no  real  observational  confirmation  of  this  fact 
exists.  Their  results  also  confirm  the  strength  of 
the  storm  although  a  microburst  impinging  on  the 
surface  would  tend  to  have  more  intense  horizontal 
outflow. 

Reflectivity  and  horizontal  wind  flow  at  4.0  km 
AGL  is  shown  in  Fig.  6.5.  This  corresponds  to  the 
level  nearest  the  base  of  the  parent  cloud.  Most 
notable  are  the  jet  axis  aligned  from  west  to  east  and 
the  zone  of  convergence  just  to  the  west  northwest  of 
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the  northernmost  reflectivity  maximum.  This  conver¬ 
gence  zone  results  from  the  flow  of  environmental  air 
around  a  barrier  formed  by  the  precipitation  core. 
The  two  areas  of  high  reflectivity  with  Z  >  40  dBz 
reflectivity  are  associated  with  the  two  main  precipi¬ 
tation  cores  within  the  storm.  The  large  area  of  20- 
30  dBz  in  the  west  northwest  corner  of  the  area  repre¬ 
sents  the  lower  edge  of  the  main  storm  cloud. 

Overall  examination  of  reflectivity  and  the  hori¬ 
zontal  wind  structure  of  this  storm  suggests  that  it 
is  in  many  ways  "typical"  of  micro  bursts,  exhibiting 
features  and  magnitudes  detected  by  other  investiga¬ 
tors  using  both  Doppler  radar  as  well  as  other  remote 
sensing  devices. 

6.1.2  Vertical  velocity 

Figure  6.5  shows  the  vertical  velocity  pattern 
and  horizontal  s t or m- r e la t i v e  wind  field  at  0.50  km. 
Contour  interval  is  2  m  s  ^  .  This  level  reveals  the 
low  level  structure  better  since  at  0.25  km  the  lower 
boundary  condition  that  w=0  at  the  surface  reduces  the 
magnitude  of  the  vertical  velocity.  Grid  distances 
are  in  km  with  the  center  of  the  grid  coincident  with 
the  microburst  at  the  surface.  Note  that  the  center 
of  the  microburst  coincides  with  the  downward  velocity 
maximum  at  this  level.  Outside  this  core  area  of 
downward  motion,  there  appears  to  be  an  annular  ring 


of  upward  motion, 
descriptions  given 


This  pattern  is  consistent  with 


in  Fujita  (19S5).  The  largest 


upward  velocities  appear  on  the  upwind  side  of  the 
microburst.  This  area  of  strong  upward  motion  results 
when  the  microburst  outflow  boundary  near  the  surface 
forcing  the  incoming  environmental  air  upward. 
Parsons  e_t  a_l_ .  (  1  9  u  5  )  also  found  upward  motion  on  the 
upwind  side  of  another  storm.  They  suggested  that  the 
outflow  from  the  nicroburst  is  negatively  buoyant  com¬ 
pared  to  the  surrounding  air.  This  dense  outward 
moving  air  would  act  in  much  the  same  manner  as  a  gust 
front,  lifting  the  more  buoyant  environmental  air 
upward. 

From  Table  5.6,  it  is  found  that  the  maximum 
average  vertical  velocity  occurs  near  the  1  km  height. 
Figure  6.7  shows  the  vertical  velocity  superimposed  on 
the  horizontal  wind  field  at  this  level.  The  annular 
ring  of  upward  motion  is  still  present  and  fairly  well 
organized.  The  horizontal  outflow  is  very  weak  at 
this  level,  again  attesting  to  the  shallow  extent  of 
the  outflow  region  from  this  microburst.  Figure  6.8 
shows  the  vertical  velocity  pattern  at  4.0  km.  This 
represents  the  vertical  velocity  near  the  cloud  base. 
With  the  exception  of  three  small  downdraft  cores 
which  are  associated  with  reflectivity  maximum  most  of 
the  air  is  moving  upward.  The  distribution  of  verti- 
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cal  velocity  at  6  km  AGL  is  presented  in  Fig.  6.5. 
Inspection  of  Figs.  6.7,  6.3  and  6.9  reveals  the 
existence  of  an  updraft-downdraft  couplet  similar  to 
that  suggested  by  Lin  e_t  a_l.  (  1985).  This  couplet, 
also  found  by  Pasken  (1981),  seems  to  be  associated 
with  more  severe  storms  and  storms  in  their  developing 
stages.  Lin  e_t  a_l .  (1985)  suggested  that  the  updraft 
is  linked  to  the  weak  echo  region  (WER)  of  inflow  into 
the  storm,  and  the  downdraft  is  associated  with  the 
precipitation  shaft  near  the  right  rear  flank  of  the 
storm.  Certainly,  examination  of  Figs.  6.5  and  6.3 
and  Figs.  6.4  and  6.6  show  that,  in  this  storm,  the 
downdraft  is  directly  associated  with  the  precipita¬ 
tion  core  of  the  storm.  Examination  of  these  reflec¬ 
tivity  and  vertical  velocity  patterns  as  well  as  other 
levels  (not  shown)  suggest  that  in  the  lower  levels 
the  strong  updraft  on  the  downwind  side  of  the  storm 
is  associated  with  lower  reflectivities  defining  the 
edge  of  the  storm  rather  than  the  WER  imbedded  within 
the  storm.  However,  the  updraft  core  to  the  upwind 
side  of  the  downdraft  is  associated  with  the  interac¬ 
tion  of  environmental  flow  and  the  storm  itself. 

The  age  and  non-severity  of  the  storm  suggest 
that  if  any  such  couplet,  i.e.,  rear  flank  downdraft  - 
WER,  had  existed  it  has  dissipated.  A  number  of 
mechanisms  could  have  contributed  to  this  dissipation. 


fwnw* 


Precipitation  loading,  as  evidenced  by  the  high 
reflectivities  at  low  levels,  could  have  destroyed  the 
updraft  region.  This  is  suggested  in  the  Byers  and 
3raham  (1949)  thunderstorm  model  where  the  updraft  is 
replaced  in  the  mature  stage  of  the  storm  by  a  down- 
draft  created  by  downward  moving  cooler  air  and 
precipitation. 

Figure  6.10  shows  the  vertical  cross  section  of 
the  wind  field  along  line  A3  from  Fig.  6.1.  The 
horizontal  component  of  the  wind  has  been  converted  to 
in-plane  wind,  i.e.  it  is  a  vector  combination  of  the 
u  and  V  component  wind  lying  parallel  to  the  plane. 
The  full  vertical  extent  of  the  microburst  is  evident. 
The  connection  between  incoming  environmental  air  at 
4.0  km  and  the  microfcurst  at  the  surface  is  clear. 

The  vertical  cross  section  of  reflectivity  super¬ 
imposed  on  the  in-plane  wind  field  is  shown  in  Fig. 
6.11.  Three  areas  of  interest  are  labeled  A,  B  and  M 
on  Fig  6.11.  Area  A  appears  to  be  a  gust  front 
associated  with  the  outflow  region  from  the  main 
reflectivity  core  (3).  The  microburst  ( M )  is  located 
in  a  region  of  high  reflectivity  gradient  and  upwind 
of  the  reflectivity  maximum.  The  horizontal  flow  is 
diverted  downward  by  this  core  of  maximum  precipita¬ 
tion  acting  as  a  barrier  to  flow.  This  effect  when 
coupled  with  downward  motion  created  by  the  precipita- 
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cion  drag  enhances  the  microburst.  Lin  e_t  a  1  .  (  1935) 
found  a  similar  mechanism  within  their  study.  Fujita 
(1985)  suggested  that  in  wet  microbursts  effects  from 
evaporative  cooling  are  reduced  and  precipitation  drag 
and  downward  deflection  of  horizontal  wind  become  more 
important.  The  updraft  on  the  upwind  side  of  the 
microburst  is  associated  with  the  edge  of  the  cloud 
boundary  distinguished  by  the  20  d3z  isoline.  It  is 
the  result  of  convergence  of  environmental  air  being 
deflected  around  the  precipitation  core  and  environ¬ 
mental  air  entering  the  storm  from  the  west  northwest. 

The  gust  front  (A)  appears  to  be  associated  with 
the  20  d3z  line  nearly  4.5  km  from  the  microburst. 
From  the  kinematic  analysis  of  Wilson  e_t  _aj^.  (  1934), 

see  their  Fig.  1,  it  appears  that  the  microburst  began 
to  form  around  1641  MOT.  The  first  sign  of  divergent 
outflow  occurred  on  the  south  side  of  the  40  d3z 
contour,  located  at  13.5  west  and  17.5  south  of  CP-2. 

Two  minutes  later,  the  divergence  increased  to 

—  3  —  i  -  £■  -  ( 

5  x  10  s  reaching  the  maximum  (2  x  10  s  )  at  1646 

MDT.  The  strong  diverging  outflow  from  the  micro burst 

center  interacts  with  the  environmental  flow  from  the 

west  to  form  a  microburst  gust  front  near  area  A  (Fig. 

6.11).  A  vortex  ring  suggested  by  Fujita  (1935)  is 

clearly  seen  between  A  and  M. 


A  snail  rotor  feature,  associated  with  the  vortex 
ring,  exists  approximately  2  km  on  the  upwind  side  of 
the  raicroburst .  This  rotor  marks  the  boundary  of  the 
outflow  region  near  the  surface  (dashed  line).  The 
vertical  effect  of  the  rotor  associated  with  this  gust 
front  appears  to  extend  to  near  the  2  km  AGL  level. 
Examination  of  other  wind  cross  sections  (not  shown) 
as  well  as  Fig  6.9  suggests  that  the  rotor  feature 
does  not  extend  completely  around  the  downburst  nor 
does  the  gust  front.  The  absence  of  the  gust  front 
a  rotor  in  certain  sections  of  the  micro  burst  region 
seems  to  be  related  to  intensity  of  the  microburst, 
motion  of  the  parent  storm,  and  the  strength  and 
direction  of  the  environmental  air  flow.  rujita 
(1935)  suggested  that  the  presence  of  a  distinguish¬ 
able  rotor  in  a  mild  microburst  tends  to  be  on  the 
upwind  side  of  the  main  outrush  center,  while  the 
feature  on  the  downwind  side  is  "smeared"  by  the 
environmental  wind. 

When  looking  at  the  vertical  velocity  patterns  in 
total,  one  gets  the  impression  that  the  storm  seems  to 
be  in  two  parts.  The  first  part  is  the  subcloud 
(below  4.0  km)  section,  while  the  upper  portion  is  the 
cloud  section.  The  prevalence  of  downward  motion  in 
the  subcloud  layer  and  only  one  portion  of  the  cloud 
section  associated  with  a  reflectivity  maximum  which 
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extends  f ron  near  the  top  of  the  cloud  to  the  surface, 
suggests  that  the  storm  may  be  separating  horizon¬ 
tally.  This  would  seem  reasonable  if  the  storm  were 
in  its  mature  or  dying  stage. 

6.3  Divergence 

Figure  6.12  depicts  the  divergence  pattern  asso¬ 
ciated  with  the  microburst  at  0.25  km.  Values  near 
the  center  of  the  microburst  are  on  the  order  of 

-2  -  i 

10  s  .  i he  area  of  convergence  w h i c h  rings  this 
central  divergence  core  is  associated  with  mass  con¬ 
vergence  between  the  outflow  of  the  microburst  and  the 
environment.  This  zone  of  convergence  is  approxi¬ 
mately  half  the  magnitude  of  the  central  divergence 
area.  In  Fig  6.12,  a  dashed  line  shows  the  microburst 
gust  front.  It  coincides  with  the  axis  of  maximum 
convergence  approximately  5  x  10~  s  Figure  6.13 

shows  the  vertical  cross  section  of  divergence  along 
line  A3  in  Fig.  6.1.  Contour  interval  is 
2.5  x  10  s'.  The  storm-relative  wind  is  superim¬ 

posed.  Divergence  exists  throughout  the  entire  region 
of  the  microburst  with  the  maximum  values  near  10  s* 
near  the  surface.  The  two  areas  of  convergence 
approximately  2-3  km  from  the  center  of  the  microburst 
represent  the  ring  of  convergence  associated  with  the 
microburst  outflow  interaction  with  the  environmental 
wind.  Maximum  values  within  this  ring  are  much 


.13  As  in  Fif>.  6.10  except  for  divergence  in 


smaller  in  magnitude  than  at  the  center  of  the  micro¬ 
burst.  Note  that  the  diverging  outflow  is  mainly  con¬ 
fined  to  the  layers  from  the  surface  to  approximately 
2  km.  Above  this  level,  convergence  dominates  due  to 
the  mass  continuity  requirement  to  maintain  strong 
downflow  within  the  downdraft  and  the  convergence  zone 
created  by  deflecting  the  environmental  flow,  see  Fig. 
6 .  S . 

6.4  Deviation  Perturbation  Pressure, 

The  retrieved  P ,  '  field  for  0.25  km  is  shown  in 
Fig.  6.14.  As  expected  there  is  a  stagnation  mesohigh 
up  to  15  Pa  inside  the  microburst.  This  stagnation 
high  is  surrounded  by  a  broad  region  of  low  pressure 
in  the  strong  outflow  area.  Such  a  pressure  pattern 
is  similar  to  those  shown  in  Fujita  (1935)  for  the 
Andrews  AF3  microburst.  The  center  of  higher  pressure 
was  located  slightly  on  the  downwind  side  of  the 
microburst.  This  is  a  manifestation  of  the  effect  of 
the  environmental  air  skewing  the  microburst  downwind. 
Fujita  (1935)  suggested  that  a  low  pressure  ring 
exists  outside  the  stagnation  high  in  a  downburst. 
This  can  be  seen  readily  in  Fig  6.14.  It  completely 
surrounds  the  storm  but  has  been  modified  by  the 
environmental  wind  flow.  The  upwind  side  of  the  storm 
has  stronger  low  pressure  t  h  a  '  the  downwind  side.  The 
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Fig.  b.14  As  in  Fig.  6.1  except  for  deviation 
perturbation  pressure.  Isobars  are  at  intervals  of 
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rotor  and  gust  front  (dashed  line)  show  up  in  the 
pressure  pattern  as  a  relative  maximum  as  one  would 
expect  from  a  convergent  situation. 

Figure  6.15  shows  that  maximum  positive  pressure 
exists  at  0.75  km  AGL  in  the  microburst.  Values  of 
+0.3  mb  are  found  in  there.  As  expected,  the  values 
are  smaller  than  those  found  in  other  studies,  e.g., 
Lin  et  al .  (1935),  Pasken  (1981)  and  Wirsing  (1985). 
However,  all  these  storms  were  much  stronger  in  nature 
than  a  microburst  producing  storm  similar  to  the  one 
being  investigated.  Parsons  ejt  2_1  .  (1935)  working 
with  another  JAWS  storm  found  retrieved  pressure 
values  of  0.5  mb  near  the  center  of  the  microburst  at 
0.25  km. 

The  horizontal  pressure  distribution  at  4  km  AGL 
is  displayed  in  Fig  6.16.  This  level  is  slightly 
above  the  cloud  base  where  the  prevailing  environ¬ 
mental  wind  begins  to  interact  with  the  main  body  of  a 
microburst-producing  storm.  Inspection  of  Fig.  6.16 
reveals  that  high  pressure,  labeled  H,  is  observed  to 
the  northwest  of  the  main  updraft  (U)  with  low  pres¬ 
sure,  labeled  L,  to  the  southeast.  The  orientation  of 
maximum  pressure  gradient  determined  in  the  vicinity 
of  U  is  approximately  in  the  direction  of  the  environ¬ 
mental  wind  shear  vector,  dV,g/dz,  from  the  northwest. 
This  shear  vector  is  determined  from  Denver's  sounding 
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shown  in  Fig.  3.3.  Our  finding  is  in  good  qualitative 
agreement  with  the  linear  theory  of  Rotunno  and  X 1 e m p 
(1982).  The  theory  predicts  that  the  dynamic  inter¬ 
action  between  the  sheared  environmental  wind  and  the 
storm's  updraft  creates  high  pressure  on  the  upshear 
flank  with  low  pressure  on  the  dovnshear  flank.  The 
magnitude  of  pressure  perturbation  near  the  updraft  is 
proportional  to  the  magnitude  of  the  height  dependent 
shear  vector,  dV^/dz,  multiplied  by  the  horizontal 
gradient  of  vertical  velocity  "FZ  w  along  the  shear 
vector.  Thus,  the  pressure  pattern  with  the  main  body 
of  the  parent  storm  displayed  in  Fig.  6.16  appears  to 
be  reliable.  A  horizontal  perturbation  pressure 
gradient,  ,  from  mesohigh  H  to  mesolow  L  is  quite 
pronounced.  The  magnitude  is  about  30  Pa  (0.3  mb) 
over  a  distance  of  5  km  which  is  much  smaller  than 
those  reported  in  studies  by  3randes  (1934),  Hane  and 
Ray  (  1985)  and  Lin  e_t  a_l.  (1985)  for  severe  convective 
storms.  It  is  worth  noting  that  the  presence  of  high 
pressure  (H)  on  the  upshear  side  of  U  plays  an  impor¬ 
tant  role  in  affecting  the  structure  and  internal 
dynamics  of  a  mic r ob ur s t- p r oduc i ng  storm.  It  will  be 
shown  later  that  this  stagnation  high  results  in 
favorable  3D  perturbation  pressure  gradients  over  that 
region  and  is  partially  responsible  for  initiating  and 
maintaining  a  downdraft  which  produced  the  microburst 


in  the  lower  layers. 

Figure  6.17  shows  the  vertical  cross  section  of 
the  perturbation  pressure  field  along  line  AB  in  Fig. 
6.1.  Contour  interval  is  10  Pa  (.1  mb).  The  storm 
relative  wind  field  is  superimposed  to  show  the 
relationship  between  the  pressure  and  flow  pattern. 
Of  interest  is  the  pressure  maximum  located  at  4.5  km. 
This  is  caused  by  the  interaction  between  the  environ¬ 
mental  flow  and  the  storm  acting  as  a  barrier.  It 
also  coincides  with  the  edge  of  the  main  storm 
deliniated  by  the  20  doz  isoline.  The  air  entering 
the  storm  is  relatively  dry.  As  it  entrains  air 
within  the  cloud,  it  becomes  cooler.  This  would  cause 
it  to  descend.  This  in  combination  with  the  inter¬ 
action  with  the  maximum  precipitation  shaft  enhances 
and  maintains  the  downdraft  which  eventually  becomes 
the  microburst.  The  ring  of  low  pressure  is  present 
and  the  low  pressure  on  the  upwind  side  coincides  with 
the  gust  front  and  rotor  near  the  surface.  The 
leading  edge  of  the  outflow  is  marked  with  a  dashed 
line . 

Fujita  (1985)  suggested  that  a  close  relation  be¬ 
tween  the  perturbation  pressure  and  the  horizontal 
wind.  Figure  6.18  shows  this  relationship  for  the 
0.25  km  level  along  line  A3  in  Fig.  6.1.  As  expected, 
and  in  agreement  with  Fujita  (1985),  when  the  pressure 
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Mg.  b.18  The  distribution  at  0.25  km  of  (a) 
deviation  perturbation  pressure  in  Pascals,  (b) 
v ertl cal  velocity  (w)  in  m  s',  (c)  divergence  in 
10  s‘o  ,  (d)  deviation  perturbation  temperature 
in  C  and  (e)  in-plane  horizontal  velocity  in  m  s* 
along  line  AB. 


is  increasing  the  horizontal  velocity  is  decreasing 
and  visa  versa.  The  asymmetry  of  the  pressure  is  due 
in  part  to  interaction  of  the  microburst  outflow  with 
the  environmental  wind.  The  horizontal  outflow  on  the 
"upwind"  side  of  the  microburst  is  slowed  down  due  to 
this  interaction  and  consequently  a  higher  pressure  is 
created  on  the  upwind  side.  This  asymmetry  is  also 
reported  by  Fujita  (19S5).  The  vertical  velocity 
profile  is  as  expected.  Upward  notion  is  associated 
with  the  convergence  and  lower  pressures,  while  the 
downward  motion  is  associated  with  the  higheer  pres¬ 
sures  and  divergence  of  the  microburst  center.  The 
asymmetrical  updrafts  are  due  to  environmental  inter¬ 
action  smearing  out  the  updraft  on  the  downwind  side 
of  the  microburst.  The  divergence  associated  with  the 
central  outflow  region  of  the  microburst  is  accom¬ 
panied  by  a  convergence  zone  on  either  side.  These 
convergent  zones  amount  to  the  micro  burst  outflow  gust 
front.  Again  the  asymmetry  in  the  central  divergence 
zone  is  associated  with  the  environmental  flow  effect. 
The  maximum  divergence  is  associated  with  the  hori- 
velocity  maximum  on  the  downwind  side  of  the  micro¬ 
burst.  The  larger  convergent  area  is  associated  with 
the  gust  front  and  the  rotor  discussed  earlier.  The 
warm  deviation  perturbation  temperature  is  skewed 
toward  the  downwind  side  of  the  microburst  center. 


155 


The  lover  teraerature  surrounding  the  central  warn  core 
is  seen  as  the  two  negative  areas. 

6.5  Deviation  perturbation  virtual  temperature, 

Figure  6.19  shows  the  deviation  perturbation  tem¬ 
perature  for  0.25  km.  Remember,  this  temperature  is 
the  deviation  from  its  horizontal  mean  instead  of  the 
deviation  from  the  environmental  temperature.  Hence, 
it  does  not  represent  the  true  perturbation  virtual 
temperature,  Ty/^’ .  It  contains  the  unknown  area  mean 
temperature  < Tv<j’  >.  For  this  reason,  caution  must  be 
exercised  when  interpreting  the  temperature  field.  A 
temperature  maximum  of  3.0  degrees  is  located  just 
downwind  from  the  microburst  center.  The  non-symmetry 
is  due  to  modification  by  the  environmental  flow. 
Fujita  (1985)  and  Srivastava  (1935)  suggested  that 
positive  deviations  of  temperature  in  this  region  are 
quite  plausible,  and  in  the  case  of  a  wet  microburst 
somewhat  more  expected. 

Srivastava  (  1985)  employed  a  simple  evaporative  1 y 
driven  downdraft  model  to  show  that  near  the  end  of 
the  life  cycle  of  the  downburst,  a  positive  tempera¬ 
ture  perturbation  developed  near  the  surface  coupled 
with  a  negative  temperature  perturbation  directly 
above,  see  his  Figs.  1  and  2.  This  pattern  can  be 
seen  in  Fig.  6.20  which  represents  the  cross  section 
of  deviation  perturbation  virtual  temperature  along 
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Fig.  6.19  As  in  rig.  6.1  except 
perturbation  virtual  temperature, 
intervals  of  1 0  C . 


deviation 


line  A3.  Srivastava  demonstrated  that  even  though  the 
downdraft  air  is  warmer  than  the  surrounding  environ¬ 
mental  air,  its  thermal  buoyancy  is  negative  since  it 
is  drier  and  hense  virtually  cooler  than  its  sur¬ 
roundings.  At  upper  levels,  the  total  buoyancy  of  the 
air  is  a  combination  of  the  thermal  buoyancy  and  the 
rainwater  mixing  ratio.  In  the  case  of  the  wet  micro¬ 
burst,  the  rainwater  loading  overcomes  the  thermal 
buoyancy  to  produce  a  negative  total  buoyancy  and 
hense  enhances  downward  motion. 

Reflectivity  patterns  suggest  that  the  microburst 
is  imbedded  within  a  precipitation  area.  Thus  air 
descending  from  the  cloud  bass  would  not  have  the 
benefit  of  evaporative  cooling  to  reduce  the  tempera¬ 
ture.  If  such  a  parcel  were  to  descend  totally  within 
a  region  of  100^  relative  humidity,  it  would  in  fact 
increase  in  temperature  in  a  nearly  dry  adiabatic 
fashion  (Srivastava,  1935).  The  second  explanation 
involves  the  age  of  the  storm.  The  value  being  com¬ 
puted  is  the  deviation  perturbation  virtual  tempera¬ 
ture.  If  the  storm  were  mature,  cooler  air  descending 
earlier  in  the  storm  could  have  cooled  the  environment 
below  the  cloud.  At  the  time  of  this  study,  the 
deviation  perturbation  temperature  would  then  be 
positive.  Fujita  (1985)  suggested  that  this  scenerio 
is  quite  possible  especially  within  wet  microbursts. 


CHAPTER  7 


Summary  and  Conclusions 

The  method  developed  by  Armijo  (1969)  was  used  to 
derive  the  3D  wind  field  from  dual-Doppler  radar  meas¬ 
urements  of  a  convective  storm  producing  a  microburst. 
The  storm  occurred  on  14  July  19S2  during  the  Joint 
Airport  Weather  Studies  (JAWS).  Vertical  velocities 
were  obtained  by  integrating  the  anelastic  continuity 
equation  downward  from  the  storm  top  near  8.5  km  to 
the  surface.  After  correcting  the  vertical  velocity 
to  conform  to  the  lower  surface  boundary  condition 
that  w=0,  the  horizontal  wind  was  corrected  using  the 
variational  approach  similar  to  that  of  Ziegler  e_t  al . 
(1983)  and  Lin  e_t  £l_.  (  1985).  This  correction  was 
necessary  to  bring  the  3D  wind  field  again  into  an¬ 
elastic  balance.  Once  the  wind  field  was  derived,  a 
thermodynamic  retrieval  technique  was  employed  to 
recover  fields  of  deviation  perturbation  pressure, 
density  and  (virtual)  temperature.  These  fields  were 
then  subjected  to  internal  consistency  checks  to  de¬ 
termine  the  adequacy  of  the  recovery  procedure  prior 
to  physical  interpretation. 
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In  an  effort  to  determine  whether  the  retrieval 
method  (which  seemed  to  perform  well  in  the  strong 
storm  environment)  would  be  sensitive  enough  to  pro¬ 
duce  usable  results  on  a  milder  storm,  a  number  of 
experiments  were  done,  see  Table  5.1.  Our  findings 
show  that: 

1)  In  order  to  extract  thermodynamic  variables  from 
3D  wind  fields,  these  wind  fields  must  be  brought  into 
anelastic  balance.  Failure  to  do  this  produces 
thermodynamic  fields  which,  at  best,  can  be  used  to 
describe  only  grossest  of  features. 

2)  Inclusion  of  a  frictional  parameterization  scheme 
in  the  computation  of  pressure  and  temperature  is 
needed,  especially  near  the  surface,  where  the  fric¬ 
tional  effect  is  most  pronounced. 

3)  Inclusion  of  a  precipitation  drag  term  is  neces¬ 
sary  for  the  temperature  recovery  scheme.  This  is 
especially  true  in  cases  of  wet  microbursts  where  this 
term  plays  a  very  important  role  in  the  retrieval. 
Failure  to  include  such  a  term  will  result  in  errors 
which  are  on  the  same  order  of  magnitude  with  the 
retrieved  temperature  and  are  biased  by  the  reflec¬ 
tivity  field. 

4)  The  thermodynamic  retrieval  scheme  is  insensitive 
to  inclusion  of  a  relatively  small  random  error  ( < 1 0% ) 
in  the  radial  wind  field.  However,  large  random 


errors  in  the  wind  field  would  distort  the  retrieved 
field  greatly,  thereby  producing  very  little  useful 
information . 

5)  Patterns  and  strengths  of  features  associated  with 
the  microburst  are  weaker  for  this  storm  when  compared 
to  those  of  more  spectacular  and  stronger  storms 
studied  by  others.  The  existence  of  a  warm  center 
below  the  microburst  was  verified.  This  conforms  well 
with  scenerios  suggested  by  Srivastava  (1985)  and 
Fujita  (1985). 

6)  Precipitation  drag  as  well  as  a  mid  level  stagna¬ 
tion  pressure  high  upwind  of  the  microburst  seem  to  be 
the  driving  force  behind  sustaining  the  microburst 
during  the  mature  and  decaying  part  of  its  life.  This 
does  not  preclude  the  possibility  that  during  the 
early  stages  of  this  event,  subcloud  evaporative 
cooling  acts  to  enhance  downward  motion  below  the 
cloud  . 

This  study  does  answer  some  of  the  questions 
raised  about  the  sensitivity  of  the  procedure  to  re¬ 
trieve  dynamic  and  thermodynamic  variables  from  dual- 
Doppler  scans.  It  also  provides  a  qualitative  if  not 
quantitative  look  at  the  structure  of  a  microburst 
producing  storm  which  was  relatively  weak  by 
comparison  to  most  of  the  other  storms  studied.  There 
are  still  a  number  of  questions  which  remain  to  be 


considered.  First,  to  what  extent  does  inclusion  of 
the  temporal  term  in  the  computations  alter  the 
retrieved  fields  and  their  magnitudes?  Secondly,  and 
from  the  aviation  point  of  view,  how  often  do  weaker 
storms  produce  microbursts  which  could  be  hazardous  to 
aircraft  performance?  Finally,  and  most  important,  is 
it  possible  to  predict  the  microburst  event  prior  to 
its  affecting  aircraft  performance  using  a  single 
Doppler  radar? 


APPENDIX  A:  Object  Analysis  Schemes 

1.  Introduction 

The  objective  map  analysis  scheme  developed  by 
Barnes  (1964)  has  been  widely  used  by  the  meteorolo¬ 
gical  community.  It  only  requires  four  passes  (3 
iterations)  to  arrive  at  an  interpolated  field  which 
fits  the  observations.  The  earlier  version  of  the 
scheme  was  modified  by  the  author  (Barnes,  1973)  to 
include  the  following  features: 

1.  It  reduces  the  number  of  passes  to  two  (one 
iteration)  to  obtain  the  desired  interpolated 
field . 

2.  It  enhances  detail  in  the  interpolated  field  by 
including  time-series  observations. 

The  first  feature  is  specifically  designed  to 
speed  up  the  rate  of  convergence,  thereby  reducing 
the  computational  time.  By  properly  changing  a 
weight  function  parameter,  details  which  require  four 
or  more  passes  through  the  data  heretofore  now  become 
discernible  with  only  two  passes. 

2.  Response  Function 

Consider  an  atmospheric  variable  f(x,y)  which  is 
uniform  in  the  Y-direction  and  is  expressed  as 

f(x,y)  ■  Asin(ax)  (Al) 

where  a(-^5^)  is  the  wave  number,  and  A  is  the 


wave 


length.  We  wish  to  obtain  the  interpolated  variable 

g(x,y)  at  some  grid  point  ( x,y )  from  the  observed 

data  at  point  (x+rcos <p  ,y  +  rsin£>  )  (see  Fig.  Al),  and 

filter  (weigh)  the  data  according  to  their  distance 

(r)  from  the  point  (x,y),  i.e., 

2TT.cc 


(  A2  ) 


where  the  weight  function  is 
v(r.k)  - 


(A3) 


and  k  is  a  parameter  (weighting  factor)  to  be 
determined . 

Equation  (3)  is  a  well  known  Gaussian  function. 
Its  shape  is  uniquely  determined  by  the  standard 
deviation  C'r  JZK  .  Hence,  larger  values  of  K  will 
give  flat  curves,  while  smaller  values  will  provide 
sharp  bell-shaped  curves. 

We  wish  to  determine  the  relationship  between 
the  observed  value,  f,  and  the  weighted  average 


value,  g,  at  the  same  point  (x,y),  i.e., 
g(x,y)  »  D(a,k)f(x,y) 


(  A4) 


where  D(a,k)  is  the  response  function  and  is  a 


function  of  the  wave  number  a  (= 


It  can  be 


shown  that  the  filtered  response  to  f(x,y)  is 

nr  ^  “CxA  -^CW3" 

D(a,k)  *  e  -  e  »  e  s  ' 


( A5 ) 


which  does  not  alter  the  phase  of  the  original 
function,  but  acts  only  to  damp  the  amplitude.  Rela- 


ri3*  A 1  Relationship  between  the  data 

and  the  grid  points  for  3arnes  scheme. 


point 
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tionship  of  the  response  function,  D(a,k),  the  wave 
length,  /X.  ,  for  various  choices  of  the  parameter  k  is 
shown  in  Fig.  A2.  It  indicates  response  correspond¬ 
ing  to  initial  (first  pass)  through  the  data  with  a 
Gaussian  type  filter  [Eq.(A3)].  In  general,  response 
is  nearly  zero  for  very  short  waves  and  approaches 
one  for  very  long  waves. 

Ideally,  one  can  select  a  small  k  which  allows 
us  to  define  the  effective  cut-off  wavelengths  (see 
Fig.  A2)  and  yet  be  confident  that  very  short  waves 
(noise)  are  suppressed.  However,  the  practical  lower 
limit  on  k  is  determined  by  the  data  density  and 
distribution.  For  example,  if  the  average  data 
spacing  is  20  km,  then  it  is  impractical  to  consider 
response  to  waves  <  40  km  in  length.  On  the  other 
hand,  if  data  are  not  uniformly  distributed,  then 
phase  changes  and  a  higher  "noise"  level  are  inherent 
to  the  interpolated  field.  Hence,  additional 
restrictions  are  required  (Barnes,  1973). 

In  his  earlier  paper  (Barnes,  1964),  the  weight 


function  hi  **  exp(-r^/4k)  was  chose  with 
'  n  i/ 


(A6) 


where  E  -  4  and  R  is  the  "radius  of  influence"  or 


simply  "scan  radius"  beyond  which  an  observation 
exerted  zero  influence  in  determining  g(x,y),  see  the 
dashed  line  in  Fig.  A3.  With  E  =  4,  we  have  repre- 


Relationship  of  response  function  to  wave 
length  A  for  various  choices  of  parameter 
K  (after  Barnes,  1973). 
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Fig.  A3  Cressnan  weighting  function  (solid 
line)  compared  to  first  pass  of  3arnes  weighting 
function  (dashed  line).  For  Cressnan,  R  *  1.7,  for 
Barnes  R  *  2.3.  Crossover  point  is  r*1.15  with 
weight  w  *  exp(-l)  *  0.353. 


sented  98 /»  of  the  influence  of  any  datum  within  the 


circular  region  whose  radius  is  R.  With  the  aid  of 
(A6),  (A5)  can  also  be  written  as 


related  to  the  mean  data  spacing  ( d^ .  It  is  also  a 
function  of  the  observational  errors  and  the  nature 
of  the  data  distribution  (Stephens  and  Stitt,  1970). 
For  the  case  of  uniform  data  distribution,  an  optimum 


choice  of  R  in  relation  to  dmis 


1  <  K/cU  <  2 


(A8) 


Figure  A4  illustrates  some  relationships  between  the 
initial  response  function  Dc  and  wavelengths  X  for 
various  values  of  scan  radius  R. 


3.  Forcing  Analysis  to  Fit  Observed  Data 

In  his  earlier  paper,  Barnes  (1964)  developed  a 
technique,  a  method  of  successive  corrections,  to  fit 
the  interpolated  function,  g(x,y),  to  the  observa¬ 
tions  within  an  arbitrarily  small  difference.  Brief¬ 
ly,  the  technique  employed  successive  applications  of 
the  same  filter  (weight  function  rj  )  to  the  residual 
differences  between  g(x,y)  and  f(x,y)  on  each  of  an 
arbitrary  number  of  passes  through  the  data.  This 
technique  was  shown  to  be  mathematically  convergent 
which  required  only  four  passes  (three  iterations)  to 
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Fig.  A4  Response  curve  for  Barnes  scheme  (Eq.  A5 
for  various  values  of  R.  Ordinate  represents  t h 
percentage  of  original  amplitude  returned  afte 
griddi.ng  procedure. 


force  the  interpolated  field  to  agree  with  observa¬ 
tions.  The  main  drawback  of  the  scheme  was  the  large 
number  of  iterations  required  to  achieve  adequate 
response  to  short  wavelengths  which  are  not  computa¬ 
tionally  economical.  Later  Barnes  (1973)  modified 
the  scheme  to  require  only  one  pass  (iteration) 
through  the  data  to  achieve  the  desired  response  at 
the  same  wavelength.  When  only  one  pass  through  the 
data  is  made,  the  initial  response  ,  is 

De  =  e  (A9) 


The  relationship  between  Dq  and  X  for  various 
values  of  4k^  is  shown  in  Fig.  A2. 

The  interpolated  field  after  the  first  pass  at 
the  grid  point  (i,j)  is  given  by 

ga(i,j)  =  f(x,y)De  (A10) 

The  subscript  d  denotes  the  first  pass  result  through 
the  data  with  weight  function  =  exp(-r  /4kft).  The 
second  pass  (first  iteration)  then  yields  smoothed 
values  of  the  residual  differences  between  f ( x , y )  and 


go(x,y)  which  are  added  to  the  first  pass 
g  0 ( i *  j )  to  give: 

g,(i.j)  =  g  Q(  i  »  j  )  +[f(x,y)  -  g  (  x  ,  y  )  ]  D, 


where 


Dj  -  exp  (-4k(  (^-)2  ) 


field  , 

(All) 

(  A 1  2  ) 


i3  the  response  resulting  from  application  of  the 
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weight  function  , 

expi-nCW) 


(  A13) 


where  k(  -  kft Y  and  0  <  V  <  1. 

Thus  t_ 

D(  «  exp((-4Yk0(%)}  =  D/  ( A 1 4 ) 

Substituting  (A10)  and  (A14)  into  (All)  yields 

g,  (  i  *  j  )  -  Dftf(x,y)  +  [  f  (  x ,  y  )  -  D0f(x,y)]DoV 

-  f(x,y)[Dtf  +  DY  -  D,DY  ] 

«  f(x,y)DJl  +  dJ"'  -  D*  ] 

-  f(x,y)D'  ( A 1 5 ) 

where  D'  =  ( 1  +  -  DY  ) 

D7-  Dd  +  (1  -  D&  ) Df  ( A 1 6 ) 

is  the  new  response  function  indicating  the  proper 

measure  of  the  degree  of  analysis  convergence.  In 

other  words,  D’  represents  how  closely  the  inter~ 

polated  values  agree  with  the  observed  ones  after  a 

second  pass  (first  iteration)  through  the  data.  The 

relationship  between  the  initial  response  function, 

D0,  and  the  final  response  function,  D' ,  for  various 

A,  is  shown  in  Fig.  A5  and  Table  Al.  Note  that  the 

curve  V  a  1.0  corresponds  to  the  response  after  one 

iteration  using  the  old  technique,  i.e., 

N  . 

DoO  »  K  Z  (1  -  D«  )  ( Al  7 ) 

The  new  technique  (kj  ■  V  kfl )  with  Y<  1  recovers 
short  wavelength  amplitudes  very  quickly.  According 
to  Barnes  (1973),  Y  cannot  be  smaller  than  0.2  with- 


t.'K 


out  creating  underflow  checks  on  most  computer 
systems.  Further,  if  'Y  >  .5,  it  does  not  produce 
rapid  analysis  convergence.  An  optimum  value  seems 
to  be 


0.2  <  V  <  0.4  (A18) 
as  discussed  in  the  study  by  Koch  e_t  a^.  (1983). 
From  (All),  the  actual  correction  pass  value  at  each 
grid  point  is  computed  as  the  sum  of  the  weighted 
averages  from  the  two  passes  with  M  observations, 
i  .  e  .  , 


I  )  tt  I  ^  ,  rrv  =  t  \r_ - - - 


S'U'j) 


2  Horn 

mu  v> 


+ 


N\ 

2 


T*\  t  \ 


Tj 


(  A 1 9 ) 


If  we  let  , 

)  ri<^=  '“r* 

then  (A19)  can  also  be  written  as 


S  *1}  ( v.y*, )  (f 1 fcu j*,) 

__  fA  -i  •  ^  Ya-=.  \  ~ — 


N\ 
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Ya 


ms  i 


o  l 
m*i 


( A19  '  ) 


where  values  of  g#(*rv?yrn)  can  be  obtained  by  a  simple 
bilinear  interpolation  between  values  of  g  (i,j)  at 
the  four  surrounding  grid  points. 


4.  Barnes'  Scheme  versus  Cressman's  Scheme 


It  is  of  interest  to  compare  the  scheme 
developed  by  Barnes  (1973)  with  the  Cressman  (1959) 


method  for  restoring  short  waves  by  successive  scans 
with  decreased  influence  radii.  Cressman's  weighting 


function , 


Wi 


is  given 

Ri-dL* 


by 


Wi  = 


r: 


>  d 


( A20) 


0  4-0 <r  Ai  >  R.fl 

where  d  is  the  distance  from  the  grid  point  (i,j)  to 
the  observation  at  (x,y)  and  RQ  is  the  scan  radius. 
The  scan  radius  Ra  (or  the  length  of  the  radius 
vector  of  influence)  is  generally  equal  to  the  hori¬ 
zontal  and/or  vertical  grid  spacing  (  )  ,  i.e., 

Rc/^  =  1  (Wilson  e_t  a_l .  ,  19S4).  Fig.  A3  illustrates 
Jressman's  weighting  function  for  R 0  *  1.7  km.  For 
comparison,  Barnes'  weighting  function,  (or  ,  is 
also  plotted  using  the  scan  radius  R  »  2.3  km 

(dashed).  Note  that  f  R.  In  this  example,  both 
curves  have  the  same  weighting  e*(»  .368)  at  the 

distance  approximately  1.15  km.  According  to  Barnes 
(1973),  his  scheme  when  compared  to  Cressman's  scheme 
has  the  following  four  advantages: 

1)  Weight  factor  4k  can  be  chosen  prior  to  analysis 
so  that  pattern  scales  supportable  by  the  data  dis¬ 
tribution  will  be  revealed  to  known  response  ampli¬ 
tudes. 


2)  Because  approaches  zero  asymptotically,  the 
influence  of  data  can  be  extended  any  distance  with- 
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out  changing  the  weight  function  and,  therefore,  the 
response  characteristics.  In  the  Cressman  technique, 
the  weighting  function  shape  is  tied  to  the  influence 
radius  (Re  )  beyond  which  zero  weight  applies.  To 
insure  that  sufficient  data  influence  the  interpola¬ 
tion  in  data  sparse  regions,  the  current  scan  radius 
is  locally  increased  in  some  applications  (Inman, 
1970)  until  a  minimum  number  of  observations  are  in¬ 
cluded.  Such  a  locally  varying  weight  factor  pro¬ 
duces  unknown  response  in  the  final  results,  and 
introduces  small  scale  irregularities  (noise)  which 
must  be  smoothed  by  later  application  of  arbitrary 
filters . 

3)  Small  scale  irregularities  are  adequately  sur- 
pressed  by  this  technique  with  the  proper  selection 
of  scan  radius  (R).  Once  the  value  of  R  is  chosen, 
there  is  no  neei  to  change  it  during  the  iteration 
process . 

4)  The  desired  pattern  resolution  can  be  achieved 
in  only  one  iteration,  instead  of  four  or  more  re¬ 
quired  with  Cressman’s  technique,  thus  effecting  a 
modest  savings  in  computer  time. 

An  objective  analysis  scheme  based  on  the  Barnes 
(1973)  technique  and  designed  for  use  on  an  inter¬ 
active  computer  was  developed  by  Koch  e_t  al .  ( 1983 )  ♦ 
In  the  paper,  a  comparison  between  the  Barnes  and 


Cressman  objective  analysis  schemes  was  discussed. 
Since  the  Cressman  weight  does  not  asymptotically 
approach  zero  with  increasing  distance  as  it  does  in 
the  Barnes  technique  (see  Fig.  A3),  it  can  cause 
serious  difficulties  when  the  data  distribution  is 
non-uniform .  Classical  sampling  theory  suggests  that 
waves  with  X<  2drn  (or  2£n)  cannot  be  resolved  pro¬ 
perly  since  five  data  points  are  required  to  describe 
the  wave  and  its  derivatives.  Both  Cressman  and 
Barnes  techniques  employ  the  method  of  successive 
corrections.  Thus,  an  adjustment  is  made  to  the 
first  pass  analysis  by  decreasing  the  scan  radius  in 
the  second  pass  through  the  data  to  recover  the 
amplitude  of  short  waves  with  wavelength  A>2drnwhich 
was  suppressed  during  the  first  pass  due  to  fil¬ 
tering.  In  the  Cressman  case,  neither  the  number  of 
additional  passes  (  4-6)  nor  the  value  of  the  scan 
radius  at  the  second  pass  are  governed  explicitly  by 
the  data  distance;  thus,  neither  is  the  filter 
response.  Stephens  and  Stitt  (1970)  developed  a 
method  of  selecting  optimum  influence  radii  for  the 
objective  analysis  of  a  scalar  field  using  the  method 
of  successive  corrections.  They  used  Cressman's 
weighting  function  to  perform  various  experiments. 
Results  showed  that  for  an  ensemble  of  randomly  dis¬ 
tributed  station  arrays,  the  optimum  influence  radius 
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occurs  at  R^/d^l.6.  However,  the  optimum  influence 
radius  increases  with  data  separation,  observational 
error,  and  wavelength  of  the  true  field  for  the 
average  taken  as  the  guess  field.  For  a  uniform  data 
distribution,  an  optimum  choice  for  the  Cressraan 
influence  radius  on  the  first  pass  can  be  made  in 
terms  of  d^.  However,  the  optimum  influence  radius 
is  not  well  known  theoretically  on  additional  passes, 
so  the  final  filter  response  is  rather  arbitrary. 

Unlike  Cress man  weighting,  the  Barnes  filter  re¬ 
sponse  characteristics  can  be  determined  prior  to  the 
analysis  so  that  pattern  scales  resolvable  by  the 
data  distribution  may  be  revealed  to  a  great  extent. 
According  to  Koch  jet  al .  (1983),  non-uniform  data 
distributions  must  be  carefully  dealt  with.  The 
average  distance  between  each  of  the  observations  to 
its  nearest  neighbor,  ,  determines  the  maximum 

detail  permitted  in  the  objective  analysis.  Once  the 
data  spacing  (An)  is  determined,  i.e.  An)^,  the 
value  should  be  properly  selected.  This  quantity 
determines  the  degree  of  convergence  of  the  inter¬ 
polated  field  towards  the  observed  field.  As  pre¬ 
viously  mentioned  (see  Fig.  A5),  the  maximum  detail 
is  obtainable  with  V*  0.2,  while  the  least  detail 
results  with  Ym  1.0.  The  small  value  of  V" 
(say  Y?:0 . 2 )  can  only  be  used  when  observations  have 


large  values  of  signal  to  noise  ratio  (SNR)  and  are 
not  substantially  contaminated  by  subgrid-scale 
atmospheric  processes.  In  addition,  the  spatial  data 
distribution  oust  be  rather  uniform.  For  the  pur¬ 
poses  of  Doppler  studies,  we  found  that  most  data 
have  a  high  enough  signal  to  noise  ratio  so  that 
can  be  chosen  between  0.2  and  0.4.  Figure  A6  pre¬ 

sents  the  relationship  between  the  response  function 
and  wavelength  for  Y =  0.3  and  R  =  2.5  km.  The  solid 
curve  represents  the  initial  response  (Do),  while  the 
dashed  curve  denotes  the  final  response  (D').  For 

y 

comparison,  the  response  curve  D0  is  also  plotted 
(dotted  curve).  It  is  seen  that  amplitudes  of  shorter 
waves  which  were  surpressed  during  the  first  pass  are 
recovered  substantially  after  one  iteration  (second 
pass).  For  example  at  A  «  3  km,  D0  *0.18  and  D*  * 
0.67  which  implies  that  for  the  wave  with  wavelength 
X-3  km,  about  67%  of  its  original  amplitude  will  be 
retained  after  only  one  correction.  On  the  other 
hand,  if  no  correction  is  made  then  only  'v*18%  of  its 
original  amplitude  will  be  kept  and  the  remaining  82% 
will  be  lost  due  to  filtering. 

From  Table  Al,  it  is  seen  that  as  the  size 
(scale)  of  waves  increases,  say  A  ^  5  km,  then  more 
than  90%  of  their  original  amplitudes  will  be  repre¬ 
sented  in  the  final  analysis  field.  It  will  be  shown 


later  that  our  final  response  curve  (D')  is  in  good 
agreement  with  the  response  curve  of  Testud  and  Chong 
(1983)  for  Analysis  B  (see  their  Fig.  2).  Inspection 
of  Fig.  A6  further  reveals  that  the  wavelengths  of 
the  minimum  resolvable  wave  for  the  scan  radius  of 
R  =  2.5  km  is  about  2  to  2.5  km  depending  on  the 
definition  of  a  cut-off  wavelength.  Testud  and  Chong 
(1983)  determined  the  cut-off  wavelength  based  on  the 
amplitude  gain  factor:  G  =  0.5,  corresponding  to  the 
response  D'  =  0.5  in  our  example.  On  the  other  hand, 
Koch  et  al.  (  1983)  used  D'  =  e"1  (=  0.378)  to 
determine  the  minimum  resolvable  wavelengths  (2  n)  in 
terms  of  known  data  spacing  An.  Consider  0.2, 
they  arrived  at  the  formula  relating  the  weighting 
factor  4k0to  the  minimum  resolvable  wavelength  2  A n , 
i  .  e  .  , 

4k0  «  5.052( ^  )“  ( A2 1 ) 
It  implies  that  the  weight  parameter  4k 0  is  fixed  by 
the  data  spacing(An)  to  give  maximum  response  of  no 
more  than  e  *  at  the  2An  wavelengths.  The  above 
relationship  will  be  slightly  changed  if  a  different 
value  of  Y  is  used.  If  we  rewrite  A(21)  as 

4k0  -  -  Cr(  (A22) 
where  Cy*  3.052  for  V  a  0.2.  Note  that  values  of  C 
become  smaller  if  larger  values  of  Y  are  used. 


It  follows  that 


BARNES  1973  RESPONSE  FUNCTION 


f.OO  8-00  12.00  16.00 

WAVELENGTH CkmD 
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8  -  n  or  =  2.87  ( A22a) 

When  ^  =*  0.3,  our  calculation  indicates  that 

C^=  3.503  and 

fc/AO  ~  2.4-  (A22b) 

which  suggests  that  the  selection  of  an  optimum  scan 
radius  is  closely  related  to  the  mean  data  spacing 
An  (or  d^ ).  If  An  =  1  km,  then  R  should  be  chosen  to 
be^-2.4  km.  On  the  other  hand,  if  the  cut-off  wave¬ 
length  is  chosen  at  G(D')  =  0.5,  then  the  relation 
(A 22b)  becomes  R/^n  ~  2.0. 

Once  the  mean  data  spacing  is  known,  the  choice 
of  grid  size  (A)  can  be  properly  made.  Since  five 
grid  points  are  required  to  represent  a  typical  wave 
on  a  grid,  and  the  minimum  resolvable  wave  is  of 
A  =  2An,  then  A  must  be  no  longer  than  one-half  of 
to  guarantee  proper  representation  of  resolvable  wave¬ 
lengths.  As  for  the  lower  limit,  the  ratio  A  /  A.  n 
should  be  less  than  1/3  to  avoid  noisy  derivative 
fields  in  velocity,  e.g.,  vorticity  and  divergence 
fields.  Koch  _et_  al .  (1983)  suggested  that 

1/3  <  1/2  ( A23 ) 

for  their  interactive  Barnes  scheme. 

5.  Objective  Analysis  Scheme  Based  on  Principles  of 
Calculus  of  Variation 

Testud  and  Chong  (1983)  proposed  a  method  based 
on  the  principles  of  numerical  variational  analysis 
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(Sasaki,  1970)  to  filter  and  interpolate  the  "raw" 
data  in  radial  velocity.  This  step  is  necessary  in 
order  to  reduce  the  noise  due  to  radar  statistical 
errors  and  to  cut  off  turbulent  motions  with  scales 
smaller  than  the  data  resolution. 

We  wish  to  find  the  filtered  function  M°  (s,l) 
which  is  the  best  fit,  in  the  least-squares  sense,  to 
the  data  set  {M^}  .  Here  s  =  x  and  1=  \]  y"2-  +  zl  in  the 
Co  plane  coordinates.  Since  M°  is  continuous  and 
differentiable  up  to  the  second  order,  the  integral 


Ton  =  r  \  it  ( tl ) : ( Ir) ' ~)  J L 


(A24) 
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is  bounded.  Here  S  denotes  the  surface  of  the  domain 

being  considered.  We  want  to  minimize  the  error 

o 

between  functions  M  (s,l)  and  Mu  such  that 

KJ 


which  satisfies  the  condition 


J,  (N  )  -  C, 

with  positive  real. 


(  A25  ) 


(A26) 


The  classical  problem  in  calculus  of  variation 
can  be  solved  by  considering  the  minimization  of  the 
function  F  in  the  form 

Fs  -jr  f  1  •+-  At  J,  Cm*  j  (A27) 
where  M,  is  a  predetermined  value.  The  Euler's 


equation  associated  with  (A27)  is 
M°  -  M  -  M.Arn  -  0 


(A28) 


v. 


-N  V.  -  A  *. 


Bn 


In  Fourier  space  (A28)  becomes 


K  = 


Kk 


(  A29 ) 


-  1-fr-UlKl 

.  .  *•  O 

where  is,  is  the  wave  vector,  and  M^,  M^denote  the 
Fourier  component  of  M°  and  M,  respectively.  Thus 
J(  (Ma  )  acts  as  a  low-pass  filter  which  has  the 
property  of  isotropy  (since  the  amplitude  gain  only 
depends  on  wave  number IKI).  Further,  the  parameter 
controls  the  "cut-off  wave  number  "1 K  |=  ]/\j  M  or  wave¬ 
length  A0-  for  which  the  amplitude  gain  (G= 

is  0.5  corresponding  to  50%  response  in  the  Barnes 
scheme  discussed  previously. 

For  the  same  reasons,  the  above  relationship  can 
also  be  applied  to  functionals  Jn using  higher  order 
derivatives ,  i . e . , 


WP. 


l 


((Zjifa(iL^)d«d<L  <A30) 
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From  Eq.(A29)  it  follows  that  the  amplitude  gain 


•  j 


G  corresponding  to  such  functionals  can  be  written  as 

=  i~mk- 

The  higher  the  order  n,  the  steeper  is  the  slope^of 
the  gain  curve  around  the  cut-off  wavenumber  IkJ—XL- 
Only  the  first  (n  ■  1)  and  second  (n  ■  2)  order 

functionals,  Jj  (M°  )  and  (M0  ),  are  considered. 
Jj (M°)  is  given  in  (24)  and  J^(M°)  has  the  form 


(A32) 


Four  numerical  experiments  were  performed  by 
Testud  and  Chong  (1983).  The  first  three  experiments 
(Analysis  A,  B,  and  C)  are  presented  below: 

Analysis  A.  -  Variational  Analysis  with  Jj. 

Consider  the  first  order  functional  Jj(M°)  with 
JAj-  0.405(km)  ,  corresponding  to  a  cut-off  wave¬ 
length  A0*=  4  km.  The  filtered  values  M^j  are  calcu¬ 
lated  from  the  matrix  equation,  see  their  Eq .(17). 
Analysis  j5  -  Variational  Analysis  with  J 

As  in  Analysis  A  except  that  Jj  (M°)  is  replaced 
by  J2(M°)  ,  the  second  order  functional.  To  assure 
the  cut-off  wavelength  Ao-  4.0  km,  pi  is  chosen  to  be 
0.164  km  . 

Analysis  C  -  Cressman’s  Filtering 

o 

Gridded  values  of  M(j  are  calculated  based  on 
Cressman's  Weight  inside  a  circle  with  a  scan  radius 
R0  ■  1.7  km.  The  value  of  chosen  gives  the  same 

"cut-off  wavelength"  4.0  km  as  in  A  and  B. 

The  filtering  characteristics  under  experiments 
A,  B,  and  C  are  presented  in  Fig.  A7.  It  represents 
the  amplitude  gain  factor  G,  defined  as  the  ratio  of 
wave  amplitude  at  output  to  wave  amplitude  at  input, 
as  a  function  of  the  wavelength.  Curves  A  and  B  are 
found  to  be  very  close  to  the  theoretical  gain  curves 


rig.  A  7  Amplitude  gain  G  of  the  filtering  as  a 
function  of  the  wavelength/..  Curve  A  variational 
analysis  with  functional  J(  ;  curve  3  variational 
analysis  with  functional  J2;  Curve  C  filterin'* 
using  the  Cressnan  weighting  function.  The 
response  curve  D'  obtained  from  Barnes  (1973) 
scheme  with  V  *0 . 3  and  R»3.6  km  is  denoted  by  the 
symbols. 
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derived  from  the  Euler  equation,  i.e., 


Table  A2  presents  values  of  G  obtained  from 


(  A33  ) 
(  A33 )  . 


These  values  are  in  good  agreement  with  curves  A  and 
B  shown  in  Fig.  A7.  It  is  apparent  that,  Analysis 
B  appears  to  be  better  than  Analyses  A  and  C.  It 
provides  a  sharper  cut-off  around  A0=*  3.6  km. 
Further,  it  completely  suppresses  subgrid  scale 
motions,  and  yet  retains  the  amplitudes  (G  >  0.95) 


for  motions  with  wavelengths  >  7  km.  Analysis  C 
appears  to  perform  better  than  Analysis  A,  but  worse 


than  Analysis  B.  For  comparison,  the  response 
function  D'  which  produces  a  cut-off  wavelength  at 

km  (  0.3  and  R«  *  4  km)  is  also  produced  in 

Fig.  A8  using  the  symbol  x.  It  is  evident  that  the 

agreement  between  the  D'  curve  and  the  B  curve  is 

remarkable . 

It  is  of  interest  to  compare  the  response  curve 
presented  in  Fig.  A6,  based  on  the  Barnes  (1973) 
scheme  with  Y»  0.3  and  2.5  km,  to  the  results 

obtained  by  Testud  and  Chong  (1983),  see  Fig.  A7.  In 
order  to  make  the  comparison  meaningful,  the  cut-off 
wavelength  A0as  revealed  in  Fig.  A6  is  now  chosen  to 
be  about  2.5  km.  Since  (33)  can  produce  the  curves 
very  comparable  to  those  presented  in  Fig.  A7  (Testud 
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and  Chong,  1983),  we  decided  to  use  (A33)  to 
calculate  values  of  amplitude  gain.  In  order  to 
produce  the  "cut-off  wavelength"  at  Xe  =  2.5  kr , 
values  of  JU,  =  0.1583  km  and  JUj_  =  0.0251  km  are 
chosen.  Results  for  Gj  and  G£  are  given  in  Table  A3. 
For  comparison,  values  for  D*  at  corresponding  wave¬ 
lengths  are  also  shown.  It  is  clear  that  G^_  is  in 
very  good  agreement  with  D'  although  the  results  were 
obtained  from  two  totally  different,  independent 
methods.  These  results  are  depicted  in  Fig.  A8  for 
illustration. 

6.  Concluding  Remarks 

The  Barnes  (1973)  objective  analysis  scheme  is 
shown  to  be  a  very  effective  method  in  reducing 
meteorological  data.  Although,  the  method  has  been 
primarily  applied  to  the  data  collected  under  the 
conventional  observational  network,  i.e.,  upper  air 
and  surface  data,  it  can  also  be  applied  to  non- 
conventional  data  sets,  e.g.,  satellite  and  Doppler 
radar  data.  The  advantages  of  the  Barnes  scheme  over 
the  Cressman  scheme  in  obtaining  the  interpolated 
field  with  successive  corrections  are  apparent.  It 
allows  the  user  to  select  a  proper  scan  radius  based 
on  the  consideration  of  the  data  distribution,  the 
observational  errors,  etc.,  prior  to  the  analysis. 
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Accordingly,  the  information  of  resolvable  waves  can 
be  preserved,  while  the  unresolvable  waves  (noise) 
will  be  properlj  filtered.  This,  it  appears  to  be  a 
useful  tool  for  Doppler  studies. 


APPENDIX  B:  Correction  of  3D  Wind  Field 


1)  Anelastic  Constraint 

The  term  "anelastic  continuity"  used  throughout 
this  study  states  essentially  that  density  is  a  func¬ 
tion  of  height  only,  or 


The  wind  field  derived  using  the  Armijo  (1969) 
equations  contain  both  real  and  erroneous  information. 
The  vertical  velocity  which  has  been  derived  by  down¬ 
ward  integration  of  the  anelastic  continuity  equation 
using  these  winds  contains  error  as  well.  This  3D 
wind  field  is,  however,  in  anelastic  balance.  The 
object  of  the  correction  scheme  is  threefold. 

1)  To  adjust  the  wind  field  so  that  the  constraint 
that  the  "true"  vertical  velocity  at  the  surface  is 


zer o  . 

2)  To  adjust  the  wind  field  so  the 
field  is  still  in  anelastic  balance. 

3)  The  adjustment  of  the  wind  field 
minimum  in  order  to  preserve  as  much 
features  as  possible. 


resultant  wind 

be  kept  to  a 
of  the  "true" 


2.  Euler-Lagrange  Equations 

Assuming  that  some  weighting  functions  G  and 
are  available  which  reflect  some  measure  of  confidence 


3. 


Proof  of  Existence  of  Euler-Lagrange  Differential 


Equation 

The  motion  of  a  particle  in  a  conservative  field 
of  forces  satisfies  both  Newton's  laws  of  motion  and 
the  Lagrangian  equations.  There  is  still  another 
relationship  which  is  satisfied  by  the  notion. 

Consider  a  given  function  F(x,y,y')  which  we 
assume  can  be  expanded  in  a  Taylor  series  in  an  inter¬ 
val  xl  <  x  <  x  2  ;  where  xl  <  x2.  The  problem  we  shall 
consider  is: 

What  function  y(x)  will  cause  the  integral  (or 


functional ) 


J=  f  F(*,y,y')  Ax 


(34) 


to  have  a  stationary  value? 

We  shall  use  the  method  of  calculus  of  varia¬ 
tions,  in  which  we  consider  neighboring  functions 
y  ( x )  +^(X)  =  F 

with  boundary  conditon 


l^(xl)  =  r^(x2)  =  0  ( B5 ) 

where  f^(x)  is  an  arbitrary  function  of  x  and  ^  is  an 
arbitrarily  small  quantity.  By  letting  ^  approach 
zero,  the  neighboring  function  can  be  made  to  approach 
the  true  function.  Figure  B.l  illustrates  this  con¬ 
cept.  The  problem  may  now  be  restated: 

What  conditions  are  necessary  to  cause  the  integral 


Fig.  31  Relationship  betveen  "true"  functional 
solution  and  neighboring  function  for  method  of 
calculus  of  variation. 


(or  functional) 


(B6) 


to  have  a  stationary  value  as  £  approaches  zero, 


i  .  e  .  , 


( B7  ) 


d>  I  ,/  =0 

The  Taylor  series  expansion  for  a  function  of  three 
variables  is 

FlV-j  ) 4- (  L +.  \ 

<)L|  ^2./' 


i  /jr  if  4  i.;'  <Tr  nri  2nF 


(  B8 ) 


where  for  n  »  2,  we  have _L_ 

Z! 


The  Taylor  series  expansion  of  (B6)  is  therefore  ob¬ 
tained  by  letting  k  »  and  d  \K'  yields 


X  (B9) 


when  (B9)  is  differentiated  with  respect  to  ^  under 
the  integral,  this  assumes  that  ^  is  not  a  function 


of  x,  the  results  are 


How  allow 


— >  0,  so  that  we  see 


(Bll) 


With  the  aid  of  the  method  of  integration  by  parts, 

& 


dx  +  liu- 


(  B 1  2 ) 


We  note  that  the  last  portion  of  (312)  reduces  to  zero 
due  to  the  boundary  conditions  we  imposed  in  (35). 

The  integral  of  (312)  is  zero  for  the  arbitrary 
function  F(x)  only  if  the  integrad  is  zero,  hence 


(313) 


which  is  called  the  Euler  -  Lagrange  differential 
equation.  The  function  y(x)  which  satisfies  (B13) 
will  cause  the  integral  of  (B4)  to  have  a  stationary 
value . 

The  resultant  Euler-Lagrange  equations  for  our 


particular  problem  are 


The  errors  in  the  observed  wind  field  cause  er¬ 


rors  in  the  divergence  field 
expressed  as 

f  (d-B) 


This  error  can  be 


(317) 


where 


is  "true"  mass  divergence 
is  observed  mass  divergen 
ji?(D  -  D)  is  mass  divergence  error. 


Integrating  (B17)  from  the  storm  top  to  the  sur- 


face  yields  / 


ZT  Zt 

But  from  the  Integral  constraint  equation  (B16) 


(318) 


one  obtains 


\N0  ~  <Z5 


(319) 


It  follows: 


( .pD<k  = 


4,' 


<2) 


(  B20 ) 

Since  the  integral  of  the  adjusted  mass  divergence 


must  vanish  we  are  left  with 

o  _  ,0  A 


J  f  J  f  0 


(321) 

'Zr 

With  the  aid  of  the  anelastic  continuity  equation 


and  the  kinematic  boundary  conditions,  it  follows  that 

a  A 


j  f  Me  =  (f  **>) 


( B22 ) 


Forcing 

the 

true  vertical  velocity  to 

vanish 

a  t 

the  storm 

top 

and  at  the  surface  is 

the 

basis 

for 

correcting 

the 

observed  wind  fields 

and 

hence 

the 

divergence  fields.  Differentiating  (B14)  with  respect 
to  x  and  (B15)  with  respect  to  y  and  add,  yields 

V_L  3X\  .  JL 


VJ-  %  J— 


( B23 ) 


In  the  above,  Dg  represents  the  vertically  integrated 
divergence  errors  accumulated  during  the  downward 
integration  process.  Since  we  do  not  know  the  verti¬ 
cal  distribution  of  Dg  well,  it  is  assumed  to  be 
independent  of  height.  This  assumption,  however, 
insures  that  the  corrected  w  does  meet  both  the  upper 
and  lower  kinematic  boundary  conditions.  Further,  the 
variationally  adjusted  three  Cartesian  wind  components 
(u,  v,  w)  satisfy  the  anelastic  continuity  equation 


everywhere  within  the  domain  being  considered.  It 
further  requires  that  the  Lagrangian  multiplier  \ 
must  be  independent  of  height.  It  follows  that  the 


analysis  equations  have  the  form: 

u-‘-  ^  7Z 

V-  V  -t- 


(324) 
( B25 )  . 


Once  the  correction  to  the  horizontal  wind  field 


is  made,  the  vertical  wind  component  can  be  recomputed 
using  the  previously  described  technique  of  downward 
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gt 4  fcW '  twx>  aw 


iffcgt  if  t  U^i.Cft.Vii 


i  *-*  •-*  *-A  l.i 


integration  with  the  corrected  horizontal  wind  field 
used  to  compute  divergence  elements. 


V.  f.  a  /.  V.  A  .•  •• 


APPENDIX  C:  Computation  of  Correction 
Weighting  Functions 


The  computation  of  the  weighting  functions  0('1~  and 
/h  used  in  the  adjustment  scheme  follows  that  sug¬ 
gested  by  Ray  e_t  a_l.  (1980)  and  Wilson  et  al .  (  1984). 
We  begin  by  applying  the  Armijo  (1969)  equations  for 
the  solution  of  horizontal  wind  to  the  standard  devia¬ 
tions  „  „ 


'^2  / 1  h  "  “ — ' 

»  \  I  _  A  v  ;  .  .  -  ' 


(ax,2.v._-  '/ 


+  c=v(u,w) 


(£*z.AV, 


+  C'W(tf.wJ  (C2). 


Wilson  et  al.  (1984)  noted  that  the  covariance 


terms  are  small  compared  to  contributions  by  the  other 
terms  and  were  therefore  neglected  in  the  computation. 
To  solve  the  above  equations,  it  was  necessary  to  add 
a  third  equation  to  close  the  problem.  A  form  of  the 
anelastic  continuity  equation  was  added  as  the  third 


equation . 

SL-  C  -  <Z 


(C3) 


This  equation  was  integrated  downward  in  exactly 
the  same  manner  used  to  solve  for  the  wind  components. 
The  boundary  condition,  *  0  at  the  storm  top,  was 


203 


assumed.  The  values  of  Ov,  andO~Vo  were  assumed  to  be 
unity  throughout  the  3-diraensional  domain.  At  first 
glance,  this  may  seem  unfounded.  The  following  test 
was  performed  to  ascertain  how  reliable  this  assump¬ 
tion  was.  At  selected  grid  locations,  the  edited 
radial  observations,  which  were  being  used  by  the 
Barnes  gridding  procedure,  were  collected  and  an  un¬ 
weighted  standard  deviation  was  computed  from  the 
radial  wind  component  from  each  radar.  The  radius  of 
influence  used  was  1.25  km.  While  means  were  under¬ 
standably  different,  these  radial  standard  deviations 
were  remarkably  well  behaved  and  clustered  in  value 
near  unity.  Recognizing  that  the  selection  of  the 
scan  radius  will  have  great  influence  upon  the  outcome 
of  this  test,  the  use  of  unity  was  not  considered 
totally  unfounded.  Wirsing  (1985)  suggested  that  the 
actual  magnitude  of  the  weighting  function  is  of  less 
importance  than  performing  the  correction  with  any 
weighting  function  and  suggested  that  perhaps  a  value 
of  unity  could  be  used  (private  communication)  and 
still  produce  credible  results. 

This  technique  has  the  appeal  of  having  computed 
the  weights  in  the  same  manner  as  computing  the  winds. 
There  are  two  primary  criticisms.  First,  it  seems  to 
some  extent  to  suggest  the  reliability  of  the 
"derived"  wind  is  more  a  function  of  position  than 


reliance  on  the  what  the  radar  "saw".  The  resultant 
field  increases  as  you  move  away  from  the  radars. 
This  would  seem  reasonable  as  the  volume  used  to 
compute  the  gate  information  increases  as  you  move 
away  from  the  radar.  Secondly,  a  correction  to 
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